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THE PROBLEM 


The National Aeronutics and Space Administration (NAS21) 
has been conducting research in digital fly-by-wire (DFBW) 
flight control for aircraft and spacecraft. At the NASA 
Dryden Flight Research Center a DFBW flight control system 
has been tested in the F-8 aircraft. In the specific flight 
control area of sensor redundancy management (SRM) a method 
capable of verifying performance through statistical analysis 
of sensor data is desirable. It is hoped that this statistical 
analysis will define performance requirements concerned with 
detecting false alarms, missed alarms, and maximum vehicle 
transients caused by sensor failures. The results should be 
used to set tolerances (thresholds or trip levels) for each 
parameter measured. 

This researcher developed a probability density function 
for the mid-value sensor select algorithm (NASA Grant NAG 4-6) . 

The probability density function was used to obtain values for 
nomograph plots for the probability of false alarm, given the 
cummulative probability of being in the domain of failure and 
system reliability. 

One of the objectives of this grant (NAG 4-9) is to apply 
the probability density function, developed in NASA grant NAG 4-6, 
to sensor output data from the AFTI F-16 aircraft. AFTI stands 
for Advanced Fighter Technology Integration and is a joint Air 
Force, NASA, and Navy program to apply futuristic aircraft 
technology using a highly modified F-16A aircraft. The AFTI F-16 
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program is expected to continue through the summer of 1984. 

More specifically, this investigator will attempt to 

. f 

reach the following goals: 

1. Rewrite the procedures, developed in NASA Grant 
NAG 4-6, for generating a probability density 
function to determine false alarm rates, using 
an algorithmic approach. 

2. Develop microcomputer software which will print 
out table of values for the cummulative probability 
of being in the domain of failure; system reliability; 
and false alarm probability, given a signal is in 

the domain of failure. 

3. Develop microcomputer software which will plot 
nomographs associated with the table of values 
in step 2. 

4. Apply the microcomputer software to sensor output 
data for various AFTI F-16 flights and sensor 
paramaters . 

5. Survey various experimental SI<M algorithms. 

6. Make practical recommendations for further research. 

Conditional probability and probability density functions 
will be used to find the probability of a false alarm, given 
the sensor select signal is in the domain of failure. Further- 
more, the resulting probabilities may be used to select 
practical threshold values. This approach also includes 
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system reliability in the determination of threshold values , 
and provides the probability of a false alarm, given an 
alarm occured. That is, if the fail light is on, the confidence 
level for resetting the system is available. The most 
difficult task will be to determine the probability density 
function for being in the domain of failure. 


One important task is the determination of the operational 
probability density function. 
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Cv=THRESHOLD VALUE 



REVIEW OF RELATED LITERATURE 


Sensor redundancy management (SRM) requires a system 
which will detect failures and reconstruct avionics accordingly. 
In the F-8 DFBW aircraft SRM is extremely important to vehicle 
performance and the activity of output signals influenced by 
device errors, failures, and vehicle motion. Failures should 
not cause excessive vehicle transient motion, which might 
result in damage or loss of the vehicle. 

Sensor redundancy management is required to detect a wide 
class of failures, detect and minimize false alarms, and detect 
whether or not the signals exceed a given trip level (threshold) . 
To maximize the detection of failures, the trip level should 
be as small as possible. However, to minimize false alarms, 
the trip level should be large enough to include the expected 
errors in the parameter channels. Errors are divided into 
three classes: fixed errors (bias errors) , noise errors 

(originating in sensor and AD converter) , time varying errors 
(dynamic errors caused by scale factor deviations, transfer 
function errors , and sensor unalignment) . 

Geiderloos and Wilson (1976) designed a Monte Carlo 
Simulation (MCS) to verify flight control SRM in the NASA 
Space Shuttle project. This MCS gives probabilities of false 
alarms, failure transients, and failure dectability as 
dicated by performance requirements. They note that MCS was 
chosen because the non-linear nature of the problem does not 
lend itself to simple linear techniques, including elementary 
statistical analysis. 


A1: the start of each run the MCS randomly chooses a number 
for each redundant parameter according to a given densicy 
function. Data is exercised through three sets of algorithms 
(Figures 1/ 2, and 3) and, for each run, finds: 

1. The maximum trip level which will give one false alarm. 

2. The maximum failure transient. 

3. The maximum trip level which will detect a given failure. 

These maxima are collected for 500 runs and then 

normalized to achieve the desired probabilities. Two of 
Gelderloos and Wilson's findings are: 

1. The probability of false alarms are a function of 
the size of errors between redundant sensors. 

2. The probability of failure detections are a function 
of the size of errors on the "good" sensors as well 
as the signal differences between the. "failed" 

and "good" sensors. 

Sensor redundancy management for the F-8 DFBW aircraft 
at the NASA Dryden Flight Research Center is divided into two 
parts (RMi and RM2) which are executed at different times in 
the computation cycle. Szalai, Felleman, Gera, and Glover 
(1976) illustrated F-8 DFBW hardware (figure 4) , software 
sequence and timing during one minor cycle = 20 m/secs 
(figure 5) , triplex SRM (figure 6) , and triplex discrete 
rendundancy management algorithm (figure 7) . A hard sensor 
fault is declared by RM2 when a sensor differs from a selected 
value by an amount greater than the allowable tolerance 
(trip level) for a given number of passes (for the F-8 DFBW 


n = given number of passes = 5=10Qm/sec) . The given number 
of passes , n, is often referred to as a window width. Figure 
8a illustrates a hardover sensor fault of sensor A, and figure 
8b illustrates a transient fault of sensor A. 

Szalai, Larson and Glover (1979) summarized the F-8 
DFBW flights in and Advisory Group for Aerospace Research and 
Development (AGARD) lecture. Table 1 lists the sensor set and 
assigned tolerances (trip levels) . Some of the findings for 
pitch rate gyro, roll rate gyro, yaw rate gyro, normal 
accelerometer, lateral accelerometer, and longitudinal 
accelerometer for sensors A and C are: 

1. RMS value of the sensor pair difference was generally 
less than half the maximum difference recorded 

(i.e. (A-C) RMS value < h maximum (A-C) ) . 

2. Maximum |a-c| < k fault threshold. 

The counters in the program logged each miscompare of 
any sensor pair, triple, or duplex set and the number of times 
the counter reached three or four (one or two of n=5) . 

Table 2 lists the cummulative sensor flight experience for 
the F-8 DFBW at NASA Dryden. 

The results are not significant because of the following 

facts : 

1. Angle of attack, pitch attitude, and altimeter 
tallies were caused by an actual hardware failure. 

2. The one lateral accelerometer fault was a false 
alarm (acceleration threshold was increased to 0.2g). 
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3. The other two lateral accelerometer and attitude 
where the count reached four were associated with 
a previous degrading of performance. 

In a Boeing Company report (1978) on Signal Monitoring 
and Voting the following reasons are given for the importance 
of sensor select and failure detection: improve auto.land 

performance, reduce nuisance disconnects, improve no failure 
performance, and add significant fault isolation capability. 

The Boeing sensor selection system is very similar to the 
Dryden F-8 DFBW system. They have also been trying to develop 
a statistical model to input sensor data, perform statistical 
analysis, identify the distributions, and predict exceedance 
rates. The Boeing Corporation has indicated that the identification 
of the distribution will help to predict low exceedance rates* 
However, they list the following as problem areas: Sensor 

data is not stationary; limited sample size (not true at Dryden 
since F-8 DFBW has had over 100 flights) ; no specific distributions 
for senior data has been identified; and, sensor distributions 
are possible to model, but extremely difficult. 

The determination of false alarm rates and practical 
threshold values would be beneficial to the AFTI F-16 research. 

Some of the futuristic capabilities of the AFTI F-16 would 
include a Voice Command System that allows the pilot and air- 
craft to converse verbally; A Digital Flight Control System 
computer and an Automated Manuvering Attack System computer 
that converse with each other to cut down pilot workload 
during combat type maneuver and direct force controls, such 


as "chin" canards that provide lift, sideforce, and axial 
translation capabilities. All of these and other features 
are designed to work together to attempt to make the 
AFTI F-16 aircraft the best fighter aircra*.u in the world. 
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METHODS AND PROCEDURES 


Given three independent signals, each measuring the same 
parameter, one of the signals or a weighted average of the 
incoming signals would be used for input in a flight control 
system. For the F-8 DFBW project, the mid-value sensor select 
algorithm was chosen: 


INPUT 1^ 

MID 


INPUT Z+ 

VALUE 

MID VALUE 

INPUT 3-j 

SELECT 



The monitoring of each of the three input signals, 
necessitates developing some criterion for determining when 
one of the input channels fails to function properly. The 
initial approach for chosing this criterion was to examine 
the difference between an input channel (say, input 3) , and 
the mid value; then assume that this difference is approximately 
normally distributed and select a value C v so that the 
probability of the "normally" distributed difference signal 
exceeding C v is quite small. 


INPUT 



MID- 


VALU| 


DIFFERENCE SIGNAL 
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As a practical matter, it was decided to determine C v 
so that the probability of n consecutive values from the 
"difference signal" exceeding C v is less than some specified 
value. The use of n values is commonly referred to as a 
window width of n. Determination of C v is predicated on the 
assumption that samples from the difference signal are in- 
dependent. The window width affects C v directly (i.e., the 
greater n, the longer it takes to eliminate a bad sensor ; 
and smaller n, the more sensitive determination of C v is to 
the kurtosis of the difference signal) . Also, the longer 
the flight time, the more likely n consecutive values from 
the difference signal will exceed C v . Thus to specify with 
at least 99% confidence, a maximum number of excursions (the 
event that n concecutive values exceed C v ) : 

i.e. p(min { |yi | }^=1— C v^ 

C v is selected so that 

p(min { |yx | >±=1—C V ) - ~ 

where T = total seconds 
n = window width 

N = number of excurisions per T seconds 
SPS = sampling rate of the system (samples per second) 

Notice, once again, that n/SPS seconds is required to 
determine whether or not the input signal is faulty. This 
time may or may not be critical. Selection of C v in this 
context simply establishes a level which the difference signal 
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is not likely to exceed (without a hard failure of some kind) . 
Figure 9 shows, for various window widths, what value C v must 
take take for having the specified number (or fewer) of 
excursions in one ifilllion frames of data. 

An alternative approach, to determining C v would be to 
consider the conditional probability of a false alarm, given 
that the difference signal has exceeded the critical value 
(C v ) . This approach is more desirable than the one outlined 
above; since it will include system reliability in the 
determination and provide the probability of a false alarm, 
given that an alarm has occured. In other words, if a fail 
light comes on, the conditional probability approach will 
provide a confidence level for resetting the system. 

The following step-by-step procedure was developed by 
the investigator. 

STEP 1 : Find P (AF/DF) , where AF=the event of a false alarm 
occuring, and DF=the event of the signal being in the domain 
of failure 



Using Baye's law we obtain: 

P (DF/AF) -P (AF) 

P (AF/DF) = — 

P (DF/AF) *P (AF) +P (DF/AF) -P(AF) 

where AF = event of no false alarm (i.e., true alarm or no alarm) 
but P (DF/AF) =1 , therefore 
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P (AP) 
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i . P (AF/DF) - — = 

P (AP) +P (DF/AF). *P (AP) 

now let A * event an alarm occurs 

AT * event a true alarm occurs 

then AF * A U AT 

and 

A A AT - + 

therefore P (DF/AF) « P(DF/(AUAT) 

* P(DFA(AUAT) ) 

P(AUAT) 

and since A/)AT * 4 

___ P (DF A) U P (DF AF) 

P ( DF/AF ) ■ 

P (A) + P (AT) 

note DF A = ♦ and DFfyAT * AT 

ljU P (DF/AF) - 

P (A) +P (AT) 

substituting ii. into i. we obtain 

P (AF/DF) = P(AF) 

P (AF) + P (AT ) P (AF) 

P (A) +P (AT) 

P (AF) CP (A) +P (AT) ] 

or iii. P (AF/DF) = 

P (AF)CP(A)+P(AT)>P(AT) -P(AF) 

Now we need formulas for P(AF) = P (false alarm, 

P(A)= p (no alarm occurs) and P(AT)=p(a true alarm occurs) 

so we can compute P (AF/DF). 
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‘‘Let P (AF) ■ P (A F) *p (alarm ftnot a hard fault) 
where A * event an alarm occurs 

P ■ event a hard fault occurs 
P * event a hard fault does not occur. 

iv. P (AF)*P(AflF)*P (A/F) •P(F)=C(DF) ' R=[/P (y) dy] • R 

DF 

where C(DF) * cummulative probability of being in the domain 
of failure 

ft = system reliability 
i.e. R » P (P) 

P(y)= operational density function 
furthermore : 

__ ORIGINAL PAGE (3 

/l»(AF)-l-P(AF) OF POOR QUALITY ’• 

P (A) = ( 1-R) +P (AP) 

v. J P (A) »1-P(A) 

V.P (AT)sl-R 


Substituting v. into iii.we obtain: 


P ( AF/DF ) * 


P(AF)CP(A) + P (AT) ] 

P(AF)[P{A)+P(AT)] + P (AT) ' P (AF ) 

P (AF) Cl~ (1-R) -P (AF) +1-R] 

P (AF) Li- (1-R) ~P (AF) +1-RJ+ (1-R) • (l-P(AF)) 

P (AF) Cl-P (AF) ] 

P (AF) [1-P (AP) ]+ (1-R) • (1-P (AF) ) 


P (AF ) 

P (AF) + (1-R) 

P (AF) 

P (AF) + P (AT) 


\ 

| 


I 


i 




l 



i 

* 


S 
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* STE*. 2: 

Determine the characteristics of the probability 
density function (pdf) * / DF P( y)dy * C(DF) 

NORMAL PROBABILITY DISTRIBUTION FOR CONTINUOUS RANDOM 
VARIABLE X(RVX): 


P(V) = 


[-1/2 (X-M) 2 /o2] 


e 


o/2rr 

where u = mean RVX 

a 2 » variance RVX 


o = st. deviation RVX 


STANDARD NORMAL PROBABILITY FUNCTION 

u = 0 a = 1 

1 -X 2 /2 

f (X) = e 

/nr 


CUMMULATIVE PROBABILITY DENSITY FUNCTION 


F (X) 




-Y 2 /2 


dy 


ORDER STATISTICS 

Now let X^, X 3 (sensor triple) 

be a random sample for RVX and 
let y^ = minimum (X^, X 3 , X 3 ) 
y 2 = mid value (X^, X 2 , X 3 ) 
y 3 = maximum (X^, X 2 , X ^ ) 

The joint probability density function of X^» X^, X^ is: 
f (X x ) -f (X 2 ) *f (X 3 ) 
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The following disjoint sets: 


*1 

= 

{<x x . 

X 2 , 

x 3 )| 

X 1 <X 2 <X 3 ) 

A 2 

= 


X 2 * 

x 3 )| 

WV 

A 3 


(<x 1# 

X 2 * 

X 3 )| 

WV 

A 4 

= 

C(X 1# 

V 

x 3 } l 

rH 

X 

V 

m 

X 

V 

CM 

X 

A 5 

= 

{(X x 

X 2 

x 3 )| 

x 3 <x 1 <x 2 } 

A 

s 

{(X , 

X , 

x )| 

X <X <X } 

6 


1 

2 

3 

3 2 1 
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are one-to-one transformations which map each of 

A^, A 2 , ... /Ag onto the same set B = (Y x , Y 2 , Y 3 ) yi^^ 

Furthermore, the inverse functions for points in: 

At are X. = y, , X 0 = y~, X, = y. 


i t 


A 2 are X^ = 


y i' x 2 = ^3 # X 3 = ^2 
A 3 are X x = y 2 , X 2 = y x , X, = y 3 


j > 

•t ' 

j * 


* 4 «e = y 2 , X 2 = y 3 , X 3 - y 3 
A 5 are x x » y 3> x 2 = y x x 3 - y 2 


it. U 3 

i. 


A areX = y , X = y , X = y 
6 13 2 *2 3 1 


The Jacobian for each of these transformations are: 


J 1 “ 


10 0 
0 10 
0 0 1 


= 1 


J 2 = 


10 0 
0 0 1 
0 10 


= -1 


i * ■ - 


f * 

5 


J 3 - 


J 5 ■ 


0 10 
10 0 
0 0 1 

0 0 1 
10 0 
0 10 


and |j^ | = 1 


= -1 


= 1 


J 4 = 


^6 = 


0 10 
0 0 1 
10 0 

0 0 1 
0 10 
10 0 


= 1 


= -1 


l = 1 / 2, 3/ 4/ 5/ 6. 
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Therefore the joint pdf of y 1# y 2# y 3 ia ; OF POOR QUALITY 

9(y r y 2 ' y 3 ) = [Jjf (y x )f (y 2 )f (y 3 ) +Jj 2 | f (y x ) f (y 2 ) f (y 3 ) 

+ *..+ J J 6 1 f (y x )f (y 2 )f (y 3 ) 
or g (y x » y 2 , y 3 ) = 6f (y^f (y 2 ) f (y 3 ) 

Now we will compute the probability density function (pdf) 
for the mid-value sensor select algorithm which is illustra- 
ted in figure 9. 

The marginal pdf for y 2 (mid-value) is : 

y 2 

h(y,) = l r 

2 ' — / y2 g(y X / y 2 < y 3 )dy 3 d yi 

y 

/-• /y 2 6f (y x )f (y 2 )f (y 3 )dy 3 dy 1 
y oo 

" «*(y 2 )/4 f (y^ /" f(y 3 )dy 3 d yi 

" 6f <v 2 ) f «y r > Ci-f(y 2 )]ay 1 

* 6f (y 2 )Cl-F(y 2 )]F(y 2 ) 
since F(X) = f f(w)dw 

THE CUMMULATIVE PDF 
and 

/.* f (w)dw = 1 

— ULI 

note that h(y 2 ) = 6f (y 2 ) Cl-F (y 2 ) ]P(y 2 ) 
has edpected value, p=0 

p (y 2 <Q)=6/2 .f ( y 2 ) Ci-f ( y 2 ) 3f (y 2 )dy 2 

since d (F (y 2 ) ) = f(y 2 )dy 2 

P(1' 2 i0 )= 6/?. F(Y 2 )df (y 2 )-CF(y 2 )] 2 df (y 2 ) 


C p (y 2 ) 3 2 _ CF(y 2 )] : 
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(0.5)2 


(0 . 5 ) 3 


or p(y 2 l°) =o*5 

Therefore our pdf has Median » 0 and is symmetric about p(mean)=0 


h(y 2 ) = 6f (y 2 )F(y 2 )Cl-F(y 2 ) 1 
h(-y 2 ) « 6f (-y 2 )F(-y 2 )Cl-F(-y 2 )D 


-X 2 / 2 


f (x) = — 


f(-y 2 ) = f(y 2 ) 

F(-y 2 ) = 1-F(y 2 ) 
l-F(-y 2 ) = F (y 2 ) 

hence h(-y 2 ) = 6f (y 2 ) [1-F (y 2 ) 3 f ( y 2 ) 

has m-G/ median = 0, and is symmetric about u. 

The variance E (y 2 2 ) -E (y 2 ) 2 =E (y 2 ) 2 
(since median = 0) 

Therefore , 

a 2 (y 2 ) variance (y 2 ) =6 J^y^ f (Y 2 ) p (y 2 ) Cl-F (y 2 ) 1 
Now we must integrate the mid -value variance. 

The formula for variance (y 2 ) above is very difficult to 
integrate . 

However, we can use the trapezoidal rule to integrate numerically 

t 

and obtain variance (mid-value)*o 2 (y 2 ) = 0.4487. 
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We define a distribution to be "near-normal if a x 2 
(Chi-squared) test of goodness fit indicates the distribution 
is indistinguishable from some normal distribution at the 
0.05 level of significance. Every normal distribution is 
"near -normal" as is h (X ) = 6f (X)F (X) (1-F (X) ) . 


The x 2 goodness fit test indicates that at the 0.05 
level of significance h(X) is not significantly different 
from a normal distribution with y=0 and o 2 = 0.^887. In 
fact, the 99 percent confidence interval for a 2 is: 

0.4222<o 2 <0.4970. 

Finally, our "near-normal" pdf for the mid-value 
selection is 

~XV2( I .4487) 

n (X) = 1 — e 

/2tF (1.2036) 

-X 2 /2 . 8974 

I J or n (X) = 0.33145 e 
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o 2 = Normal Part + Near-Normal Part = 1+0.4487=1.4487 


a = /l. 4487 = 1.2036, and 

-3o <_ 99 percent of signals <_ 3a 
-3.099 <_ 99 percent of signals <_ 3.099 

STEP 3 ; The present investigator modified the NASA Dryden 
prepackaged programs FLIFRNT (for obtaining sensor channel data) 
and SPA (for obtaining a statistical analysis of sensor -channel 
data) to obtain this information for differences (left-self, 
left-right, self-right) between sensor channels (see Appendix C) . 

STEP 4 ; Original microcomputer software (see Appendix D) was 
written to obtain tables of values and plots of corresponding 
nomographs for P (AF/DF) , given the sensor channel difference 
was in the domain of failure and system reliability. 

STEP 5 ; The revised SPA program and original microcomputer 
software was run using sensor data collected from various 
AFTI F-16 flights. 


FINDINGS 


, t 
A <-• 




The NASA Dryden computer center has a statistical 
package, entitled SPA (see Appendix E) , which can cause linear 
trends to be removed from unfiltered data; filter data; 
perform descriptive statistical analysis; nonparametric 
statistics, root mean square analysis, spectrum analysis, 
histogram plots, and normal curve fitting to histogram plots. 
Figure 10 (roll rate gyro sensor data for channel C on flight 
22) and figure 11 (longitudinal axis sensor data for channel 
C on flight 22) are indicative of F-8DFBW sensor data dis- 
tributions. Most sensor data distributions for the F-8DFBW 
aircraft have the following properties: the data is subject 

to abrupt peaking; few or no extremity values exist; and observed 
chi-square values exceed critical values (see Table 3) . These 
non-normal properties induced the investigator to formulate 
the near-normal probability density functions in the previous 
section of this paper. 

Subsequently , for NASA Grant NAG 4-6 a BASIC computer 
program was constructed and executed using sensor data from 
roll rate gyros A and C for flight 23 of the F-8DFBW aircraft. 
This program accepted a sensor value; then computed the 
corresponding value of C(DF), using the probability density 
function; then let R = reliability of the system assume values 
from 0.9999 through 0.9950; and finally compute the corresponding 
values of P (AF/DF) probability of a false alarm, given the 
sensor select value was in the domain of failure, 
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The results of this BASIC program were used to construct 
the general nomograph, illustrated in figure 12. As an 
illustrative example, we see in figure 12 that for 
C (DP) = 5 . 2xl0 -5 and R = 0.999, P(AF/DF) = 0,05. 

Now, introducing the selection of a desired window 
width (which was explained in the previous section of this 
paper) ; to the BASIC program and general nomograph which 
computes P (AF/DF) , given C(DF) and R, a general nomograph 
for selecting C v (critical threshold values) was constructed 
in Figure 13. For example, figure 13 illustrates that for a 
window width of n = 5, the threshold value 
C v s 1.48a s= 1.48(1.2036) s 1.78 (for the mean-value sensor 
select algorithm) . 

Repeating the processes outlined above the investigator 
obtained a second output of values for the mid-value sensor- 
select pdf. The statistical comparison of the mid-value and 
mean-value sensor select distribution, as expected, was not 
significant at the 0.05 level. The general nomographs in 
figures 12 and 13 apply to both mid-value and mean-value 
sensor select probability density functions. 

Modified versions of the FLIFRNT and SPA computer 
programs (Appendix C) , and microcomputer software (Appendix D) 
were used to obtain complete statistical analyses, false 
alarm probabilities and corresponding nomographs for AFTI F-16 
flight data collected by the investigator, during August, 1982 
(Appendix F) . 
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A COMPARISON OP VARIOUS SRM ALGORITHMS 


After one sensor has failed in a triplex system, the 
SRM algorithm must perform selection, detection, and recon- 
figuration on the remaining two sensor inputs. 

One basic idea in SRM, for a triplex system, is that 
when a particular sensor input differs from the other two; 
it is probably true that this particular sensor, and not the 
other two, has failed. 

A miscompare between two inputs creates a problem for 
the computer, because it has no way of telling which input is 
correct. The computer must have a third source of information, 
before it can isolate the bad input and avoid declaring both 
inputs bad. 

(A) ASPECTS OF SRM 

1. Selection Process 

2. Fault Detection Process 

3. Failure Reconfiguration Process 

(B) SELECTION METHODS FOR SRM 

1. Mid-Value Select (MVL) - takes the three 
input signals and picks the value "inbetween” 
the other two (see figure 14) . 

2. Self-If-Good (SIG) - each processor uses its 
respective sensor input, providing it was not 
previously declared out of tolerance (see figure 15) . 

3. Averaging (AVE) - mathematical average of the 

three inputs (see figure 16). Notes: if one sensor 
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has failed , the remaining two are averaged; and 
only the MVL method allows selection to be run 
before detection. 

(C) DETECTION METHODS FOR SRM 

1. Result Minus Sensor (RES-X) - takes the result 
of the selection process and subtracts each sensor 
value for comparison to a tolerance (see figure 17) . 

RES-X is performed every 20 m/secs and each miscompare 
represents one persistence count (PC) . Given a 
frame rate and a set time, a fault must persist for 

a specified maximum PC before being declared permantently 
failed. 

2. Differening (DIPP) - compares the difference of 
each sensor pair to a tolerance for determination of 
a miscompare. If two comparisons fail, the sensor 
common to those tested has its PC incremented. 

NOTES : 

a) Acceptable selection/detection processes: MVL/RES-X, 
MVL/DIFF, SIG/DIFF, AVE/DIFF. 

b) RES-X requires selection first, and only MVL allows 
selection before detection. 

(D) RECONFIGURATION METHODS FOR SRM 
Reconfiguration indicates what must be done when only 

two good sensors remain and a miscompare between these two 
occurs. The needs for reconfiguration are to set system 
reliability requirements and indicate the importance of the 
sensor in question. 


1 




. Resonability Check - determines if the input is 
within its physical range. For example, if an air- 
craft cannot reach an altitude greater than 4 OK, a 
value greater than 4 OK would be used as a resonability 
check for altitude. 

2. Rate Change - tests the sensor input in question 
against a known physical limit, namely a given rate 
change. The rate of change and resonability methods 
are quite similar. 

3. Analytic Redundancy Management (ARM) -makes use of 
unfailed sensors which are related to the sensor being 
monitored. The ARM relationship is modeled in the 
digital system and driven by dissimilar sensor inputs. 
The main function of the ARM algorithm is to output 

a third sensor value, in order to isolate the fault. 

4. Failure Analysis - uses knowledge of the system hard- 
ware structure to isolate failures. If a number of 
analogue inputs fail at the same time; a failure 

of analogue to digital is declared, then all analogue 
inputs for the channel in question are put in a 
failed status. If a second miscompare occurs and 
this failure is analyzed as a failure of a higher 
order device (such as analogue to digital converter) , 
then the corresponding input would be declared failed, 
allowing the remaining good sensor to be identified. 


RECOMMENDATIONS 


It seems that conditional probability and order statistics 
offer an elementary, interesting, and promising approach to the 
complex problem of analyzing sensor distributions, predicting 
false alarms, and selecting practical threshold values. 

In conclusion, the investigator offers the following 
specific recommendations for further research: 

1. Use the probability function and microcomputer 
software developed in this research paper to perform 
a statistical analysis for each of the sensor 
parameters on current APTI F-16 flights. 

2. Attempt to use the statistical analysis in step 1 
to set practical threshold values for APTI F-16 
sensor parameters. 

3. Develop microcomputer software to replace the NASA 
Dryden SPA statistical package. 

4. Develop statistical procedures to study the various 
experimental SRM algorithms. 

5. Develop statistical procedures to compare the results 
of various experimental SRM algorithms. 

6* Determine the probability of a sensor miscompare 

based on noise frequency content and sampling delays 

in an asynchonous system. It is possible that an 

asynchonous system will reduce the input of electrical 

transients on system operations. One theory is that 

if an electrical transient occurs, it will .appear on 

only one channel at a time and will be voted out in 
the selection process. 
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(a) Major hardware elements. 
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(b) Pallet assembly . 

Figure 4. BFBW hardware elements with detail of pallet assembly . 
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Selected value 


Figure 14. Analog eeneor fault detection isolation and reconfiguration algorithm. 
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Figure 16 , Mean-Value Sensor Select Algorithm 
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Figure 17. RES-X Detection Method for SRM 
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SENSOR SET 

TOLERANCE 

Pitch rate (deg/s) 

4 deg/sec 

Roll rate (deg/s) 

10 deg/sec 

Yaw rate (deg/s) 

4 deg/sec 

Axial accelerometer (g) 

O.lg 

Lateral accelemoter (g) 

0.2g 

Normal accelerometer (g) 

0. 5g 

Pitch C stick (cm) 

1 . 0cm 

Roll -C stick (cm) 

1.0cm 

Rudder pedals (cm) 

0.5cm 

Angle of attack (deg) 

2 . Odeg 

Left secondary actuator (cm) 

3 . 5cm 

Right secondary actuator (cm) 

3 . 5cm 

Pitch attitude (deg) 

15deg 

Roll attitude (deg) 

15deg 

Heading (deg) 

15deg 

Mach number 

0.05 


Altitude (m) 


L50m 




SENSOR 

NUMBER OF 
ACTUAL 
HARDWARE 
FAULTS 


YAW RATE GYRO 

0 

YAW RATE GYRO 

0 

YAW RATE 

0 

LONGITUDINAL ACCELEROMETER 

0 

LATERAL ACCELEROMETER 

0 

NORMAL ACCELEROMETER 

0 

ALTIMETER 

1 

MACH 

0 

PITCH ATTITUDE 

1 

ROLL ATTITUDE 

0 

HEADING 

0 

PITCH CENTER STICK 

0 

ROLL CENTER STICK 

0 

RUDDER PEDALS 

0 

PITCH SIDESTICK 

0 

ROLL SIDESTICK 

0 

ANGLE OF ATTACK 

0 


Table 2: F8 DFBW Sensor Flight Experience 


NUMBER CF 

FALSE 

FAULT 

DECLARATIONS 

0 

0 

0 

0 


NUMBER 
OF TIMES 
FAIICOUNT 
REACHED 4 


0 

0 

0 

0 


MAXIMUM 
NUMBER OF 
MISCOMPARES 
ANY FLIGHT 


0 

0 

0 

0 


1 


226 



6 

267 

0 


0 


2 


0 


0 


4 


' 0 


0 


0 


10 

72 

9 

0 


0 

0 

0 

0 


0 

0 

0 

0 


0 

0 

0 

0 


0 


0 


7 













CRITICAL VALUE 


X 2 VALUES 

_ - x 2 - 

df 

(0.995 level) 

QBA-QBC 

16475. 

8 

21.96 

PBA-PBC 

186386. 

10 

25.19 

RBA-RBC 

117647. 

14 

31.32 

NXA-NXC 

200559. 

24 

45.56 

NYA-NYC 

465560. 

27 

49.64 

NBA-NBC 

491741. 

21 

41.40 


NOTE: x 2 VALUES COMPUTED UP TO EXPECTED FREQUENCY 

LESS THAN 5. 


IN EACH CASE COMPUTED VALUE EXCEEDS CRITICAL VALUE SO 
DISTRIBUTIONS MAY BE TAKEN' TO BE DIFFERENT FROM NORMAL. 


Table 3: Chi-Square Test for Sensor Select Values 

(Mid-Value Algorithm) 
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APPENDIX 
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MODIFIED FLIFRNT PROGRAM 
MODIFIED CARD SETUP FOR SPA PACKAGE' 
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1*1 

102 

10 • 
2 u . 


0 <1* 101) TITLE 
FURMAT (8A10T 
IF(E0F< 1 ) 

WRITE <2) 

READ 102» UFLID 
FORMAT ' R 1 0 ) 

RERD FLIN 
PRINT FLIN * 

PRINT 103» UFLID 
FORMAT (IHO. "UFLID 


U IT f r *3*c »;•** • • 


R10T — FWI - »i '»'■ * ■ ' , » p p i ’! .1 ■ ymn - 

.NE.0> CO TO 5 1 \ *' * 1 ||* 

TITLE • . • - X \W L> - 

|iP| >f| . 7»> -A . ^ UVA.Uy~'fc« V , . fc-L. .... .>■»*• . • 




202 

203 

2 *' 


105 

2 


J " 
\ 


ORIGINAL PAGE LI 
, 0F p 00R QUALITY 

• RIO. ') 

IF(NCH.LT.l .OR. NCH.GT.10) PRINT 201 ' \ • ! 

FORMAT < IX. "REQUESTED NCH IS OUTSIDE OF RLLOWRBLE LIMITS. CORRECT 
-AMD RESUBMIT. " ) 

I F (NCH. LT . 1 .OR. NCH.GT.10) GO TO 999 
WRITE <2> NCH.SPS 

"ERD 101* (PARMI D< I ) » I * 1 » NCH) * * 

° I TE <2) <PARMID< I). I=1»NCH) .• 

. NOT. PR I NT I T> PRINT 104. (PARMI D< I ) . I » 1 . NCH > 
r RMAT(1X." PR PM IDS REQUESTED RRE: ".iOAlO./- 
ISHR = ST(1) » ; 

ISMIN = ST<2) . . 

I SEC * ST< 3) 

ISMIL » ST( 4 ) 

IEHR « ET(1) 

I EMIN * ET<2) ' 

IESEC = ET<3) , . - 

I EMIL = ET<4) . 

CRLL DRBOP < IUCD. UFLID. NCH. PRRMID. MEAS. IEU.CALID. ICV. PRRMNRM) 

IF* ICCPRM. NE. 0 > PRINT 202 

FOPMRT <. 1H0.-RT LERST ONE REQUESTED PARAMETER IS MISSING FROM THE 
-FLIDRB. CORRECT AND RES'JBMIT, “ ) 

I F ( ICCTME . NE. 0 ) PRINT 203 ‘ *V 

FORMRT < 1 HO. "AT LEAST PART OF THE . REQUESTED TIME IS MISSING FROM T 
-HE FLIDAB. CORRECT AND RESUBMIT.,") « .* > 

-( I CATS. NE. 0) PRINT 204 ' - - 

ORMAT < 1H0. "CATASTROPHIC FLIDRB HEADER ERROR.”) 
i F < ICCPRM. NE. 0 .OR. ICCTME. NE.0 *0R, ICATS.NE.O) GO T0.9SS.- 
RST<1) = ISHR " ' 

fiST*2> = ISMIN . ~ ✓mf.V ** 

AST (3) * I SEC • •> 

AST <4 > = ISMIL * . : # , r . r 

IF'. PRINT TT) PRINT 105. TITLE. (PARMID( I ). >1 . NCH) 

FORMAT < 1H1 » 8R1 0. /VT15. "UNFILTERED DATA FROM , FLIDAB" » -VT17. 

- 9(A?« 5X) . R7) . • r * 

CALL DABIR (MEAS.DAT, ITYPE. ICV. NCH) 

IF( ICCDRT . GT. 0 ) GO TO 3 • 

MILL I = TT.TT. 

WRITE <2) I D ITT D2 » M I LL I » < DRT < I ) » I * l » NCH ) J \ 


t I < » 

I i* 

" * 


1 0<i 


■pATTI * 

j>* a. * •5>*rfo 

PATS * 

! Dfnfa * Vfrr\ 
>rrh) * 


U| 


999 


I F( PRINT IT) PRINT lOSr *ID1 . ID2. (DflTCI ) , >T»*NCtf>*' 

FORMAT < IX. A 10. A2. 10G12. 5) 

GO TO 2 . . i. 

ENDFILE 2 . 

AET ( l '* = IDHP jf a ^ ■ 

A" '<2 ' - I dm I n — ** ^ »— — $rf ^ 

A.-.U3) = IDSEC / • $=. (p^d-^3 

f iGT <4 > = IDMIL / $ r 

1-^ I NT 107. AST. AET 

FORMAT, (1H0. "ACTUAL START TIME = f* . 12. IX. 12. IX. 12. IX. 1 3. 10X. 

* "ACTUAL END TIME = ". 12. IX. 12. IX. 12. IX. 13) 

• RITE (2) AST .AET i 

GO TO 1 
ENDFILE 2 

REWIND 2 1 ! , * 

STOP 
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APPENDIX D 

i 

MICROCOMPUTER SOFTWARE FOR COMPUTING j 

FALSE ALARM RATES AND PLOTTING 
CORRESPONDING NOMOGRAPHS 



10 

20 

30 

40 

$0 

4 

71 

80 

85 

90 

100 

110 

145 

150 

160 

170 

180 

190 

220 

225 

230 

325 

330 

328 

345 

350 

388 

501 

502 

504 

505 
510 
515 
320 
525 
530 

3 

- ,0 
545 
550 
555 
560 
562 

565 

566 

567 

568 

570 

571 

30 
585 
590 
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REM 

REM — NASA PROGRAM — WRITTEN BY JEhrfY H. SITBON 
REM -- PROGRAM WRITTEN ON NORTH*STAR COMPUTER <HP PAS l 

REM 

REM 

REM SELECTION PROGRAM 

DIM S9*<255) ,F2*<8) ,F*<33> .FI *(33) , A*( 15) 

DIM 06 ( 1 0000 ) >P(2t 300 ) , Z9 * < 10) 

GOSUB 10000 \ REM LEARN TERMINAL CODES 
PRINT -'I* 

GOSUB 9500 \ REM CHRCK FOR DIR FILE AND DEFAULT DRIVE 
PRINT Cl*. Bl*. HI* \ REM CLEAR SCREEN, RING BELL, HIGH 
PRINT TABI20) , "SENSOR REDUNDANCY MANAGEMENT" 

TAB ( 20 ) ,FNS*<28,45> 


T* ( 33 ) ,0(3,20) , P* ( 25 ) 


INTENSITY 


PRINT 

Z»5 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


\ PRINT 
TAB ( Z ) , " 
TAB ( Z ) , " 
TAB ( Z ) , “ 
TAB ( Z ) , " 
TAB ( Z ) , " 
TAB ( Z ) , " 
TAB ( Z ) , " 
\ PRINT 
" T v r e in 


PRINT 
, H 1 * , 
, H 1 * , 
,H1 * , 
, H 1 * , 
» H 1 * , 
, H 1 * , 
.HI*, 


PRINT Nl* 


N 1 * i 

” 2 “ , N 1 * , 
“ 3 " , N 1 * . 
" 4 " , N 1 * , 
"5" ,N1*. 
" 6 " , N 1 ♦ , 
" 7 " , N 1 * . 


, TAB ( 15) 
TAB (15) 
TAB < 15) 
TAP (15) 
TAB( 15) 
TAB ( 15) 
TAB( 15) 


"Enter data to Pile" 

"Display data from Pile" 

"Edit data in Pilt" 

"Compute and ararh probabilities 
"Delete Pile" 

"Display Pile list" 

"Exit prosram" 


the number oP your choice 


mm A mM m c m 335 
PRINT A* 

ON A-4B GOTO 500,900,90,3000.90,6000,8990 

B 6 S°e 89 er DATA TO FILE 

PRINT Cl*. .1* ,B1*,CHR*( 10) ,CHR*( 10) ,CHR*( 10) 

T A r .5), "INPUT MODE" \FOR R*1 TO 1000 \ NEXT R 
50o \ GOTO 575 
Cl*. Bl*. HI* 


PRINT 

GOSUB 

PRINT 

2=15 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


TAB ( Z) , "SENSOR REDUNDANCY MANAGEMENT" 
TAB(Z) ,FNS*<28»45) \ Z=ll 
PRINT 


TAB(Z-*-5 ) , "Probal i I i t y oP Palse alarm" 
PRINT TAB ( Z+5 ) ."Enter data prosram" 

PRINT TAB(Z+5> , "Written by J. Sitbon" 

T AB ( Z+5 ) , “For F. Scalzo Ph.D" 


TAB ( 8 ) , " TYPE IN THE NASA PROJECT FILE NAME" 
\ PRINT TAB ( 10 > , FNS* (25,45) ,U1*,U14,H1* 


project Pile r.ame' 
(Yes or No) ", 


DIRECTORY" ,B1* 


• Bl* 


605 
610 
615 
620 
621 
622 
62 I 
625 
630 
635 

113 

651 

660 

665 

0 

5 

680 

700 

705 

707 

708 
7 09 
710 


\ 

\ 

\ 

\ PRINT 
Nl* 

\ PRINT 
\ PRINT 
TAB ( 10) , 

F* ( 1 , 33 ) =FNS* (33,32) 

INPUT 1 " " , F* ( 1 ,25) 

A**FNT* ( F* ) 

IF L=0 THEN GOTO 565 

PRINT \ PR I NT \ PRINT Nl*\ PRINT TAe ( 1 0 ) , "LOOK I NG UP DISK 
GOSUB 2000 \ REM SEARCH DIR 
RETURN 

^RINf 1 \ T P^IN?°Cl*'f^NT*(F*) , " is a new 
PRINT \ PRINT 

PRINT "Should this Pile be created 7 
INPUT 1 " " , A * 

IF LEN ( A* ) =0 THEN GOTO 595 

IF A * ( 1 , 1 ) = " N " THEN GOTO 90 \ REM GOTO MENU 

IF A*( 1 , 1 )<>"Y" THEN GOTO 580 

G0SU3 2500 \ REM WRITE DIR ENTRY 

PRINT \ PRINT 

IK M*0 THEN GOTO 625 

FOR R= 1 TO 500 \ NEXT R 

GOTO 90 

PRINT "Project Pile ",T*»" will be stored in target file ",F2*(1,8) 
PRINT 

PRINT T AB (10) , "PRESS ANY KEY TO CONTINUE ",B1*, 

AS| N jNCHAR*(0) \ GOTO 700 

PRINT Cl * , Bl * , " Pro Jec t name ",T*>" exist, 

PRINT \ PRINT "SHOULD THIS FILE STORE NEW 
INPUT 1 " " , A * 

IF LEN ( A* ) *0 THEN GOTO 665 
IF A* ( 1 , 1 ) = "N " THEN GOTO 90 
IF A* ( 1 , 1 )< >"Y" THEN GOTO 650 
REM INPUT DATA 
LINE#0,80, 0 
PRINT Cl*. HI *,31*. 

X=LEN ( T* ) + 1 S 
A = I NT ( 36-X/ 2 ) 

PRINT FNP* ( A , 0 ) , T* , " (INPUT MODE)’ 


stored 
DATA " 


in target Pile ",F2*(1,8) 
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( INPUT MODE) ' 


713 

720 

723 

730 

731 
735 
VO 

5tiO 

770 

780 

790 

800 

810 

820 

133 

850 

8G0 

B* 

901 
903 
90G 
910 
920 
930 
94 0 
950 
960 
970 
980 
905 
930 
: 

18*3 

1030 
1040 
1050 
JO GO 
70 
1-..75 
1 0 -''to 
1080 
1085 
10, 

1 1<: 
1110 
1120 
1 1 30 
1 300 
1310 

1330 

1331 

1332 
1350 
1353 
1 3S0 
1370 
1380 
1385 

1400 

1410 

1411 
14 20 
14 30 
14 40 
1450 
2000 
2003 
2010 

|«1? 
1 *20 
2025 
2030 
2035 
2040 
2500 

m 


PRINT 

PRINT 

PRINT 

gosub 

PRINT 

gosub 

PRINT 


FNP* ( A « 1 ) » FNS* ( X » 45 ) 


TAe< 13) 
FNS* <75 
1010C \ 
HI*. 
10500 \ 
C1*.N1* 
LINE #0.80 
PRINT "SAYING 
PRINT \ PRINT 
OPE* #1. 

FOR Y ■ 1 
FOR X» 1 
WRITE #1 


•L-S" . T AB ( 38 ) , 
45) ,N1* 

REN DRAW GRID 


•L-R M .TAB(63) ,“S-R' 


REM 

Bl* 


INPUT DATA TO GRID 


ORIGINAL PAGE 13 
OF POOR QUALITY 


STAND BYE' 


#1 

\ 

'1 


PRINT "DATA SAVED 
TO 300 \ NEXT R 


.61* 


FILE " » T» » " ON DISK" 

\ PRINT TAB ( 25 ). "PLEASE 
F2*< l , 8 ) ♦ " . " +STR* ( D ) 

TO Y6-3 
TO 3 
, D(X.Y) 

NEXT 
NEXT 
CLOSE 
PRINT 
FOR R< 

$?M°R?ADING DATA FROM FILE 
PRINT C1*.H1*,B1*,CHR*< 10) .CHR*<10> .CHR*< 10) 

PRINT TAB< 15) . "#*#*** DISPLAY DATA FILE #*•*****" 

FOR R* 1 TO 1000 \ NEXT R 
GOSUB 503 \ REM GET FILE 
IF M * 1 THEN GOTO 985 

PRINT \ PRINT TAB(20) . "Can 't find " , FNT* < F* < 1 , 25 > ) \ PRINT 
PRINT "Would you like to try aaain? ". 

INPUT 1 " " . A* \ IF LEN(A*>»0 THEN GOTO 950 
IF A * ( 1 » 1 ) * " N " THEN GOTO 90 
IF A* ( 1 » 1 ) * " Y" THEN GOTO 910 
GOTO 930 

GOSUB 990 \ GOTO 1080 

PRINT \ PRINT "FOUND THE FILE “ » T * » " . READING IN DATA" 
OPEN #1 ,F2*< 1 ,8) + " . "+STR*(D) 

FOR Y-l TO ;0 
FOR X»1 TO 3 
READ# 1 » D(X.Y) 

NEXT 

IF TYP ( 1 ) »0 THEN EXIT 1075 

NEXT 

Y«Y-1 

CLOSE# 1 

RETURN 

S«0 \ REM DISPLAY 

GOSUB 1300 \ REM DO DISPLAY 

PRINT FNP* ( 22 . 23 ) » "H = HARDCOPY E * EXIT ", 

A* » I NCHAR* ( 0 ) 

IF A * » " E " THEN GOTO 90 
IF A * < > " H " THEN GOTO 1100 
S*1 \ GOTO 1085 
REM DISPLAY DATA ONLY 
IF S»0 THEN GOTO 1350 
REM PRINTER 

FILL 65534,84 \ FILL 65532, 88\S=1 
•#S CHR*<27) , "COw" ,CHR*<27) ,"C2z" 

A = INT (36-LEN(T*)/2) 

IF S= 1 THEN PR I NT# 1 CHR* ( 10 ) » CHR* ( 1 0 ) , CHR* (10) 

IF S»1 THEN PRINT #1,TAB(A),T* 

IF S*0 THEN PRINT #0 H 1 * . Bl * , Cl * . FNP* ( A , 0 ) . T* . Nl * 

PR I NT #S TAB ( A ) r \A«LEN ( T* ) \ PRINT#S FNS*iA,45> 

IF S*1 THEN PR I NT # 1 CHR* ( 1 0 ) . CHR* ( 1 0 ) 

PR I NT#S TAB( 13) . "L-S" , TAe<3B) , "L-R" , TAB ( 63 ) . "S-R" 

IF S*1 THEN PR I NT# 1 FNS*?75 . 45 ) . CHR* ( 10) 

FOR X* 1 TO Y 

PRINT #S X,TAB(7).D<lfX),TAB(32)»D(2*X)/TAB(57).D(3,X) 

IF S*1 THEN PR I NT # 1 
NEXT 

IF S*0 THEN RETURN 
A*CALL ( G5 1 G9 ) 

RETURN 

REM SEARCH PROGRAM OF DIR FILE 
OPEN #17.5, "DIR, "+STR*(D) 

M*0 \ P=-25 \ REM NO MATCH 
<1)=0 THEN GOTO 2035 


KpIJS* 


READ # 1 . F 1 * ( 1 . 33 ) 

IF Fl*( 1 ,25) =F*( 1 .25) THEN M*1 ELSE GOTO 2015 
F2*( 1 » 8 > *F 1 * ( 2G , 33 ) 

CLOSE #1 
RETURN 

REM WRITE PROGRAM TO DIRECTORY 
M*Q \ REM DISK NOT FULL 
A = G5 


- 52 - 


2508 F2*( l ,8) *FNS*(8.32> 

2510 F2* ( 1 » 8 ) ■ "NASA#-- " +CHR* ( A ) 

2511 A*»F2*+","+STR*(D) 

2515 IF FILE(A*)«-l THEN GOTO 2T30 
2520 A* A* 1 
2.525 GOTO 2510 

30 OPEN #17.5. “DIR. " + STR*(D> 

1.35 P*0 \ REM POINTER 

2540 IF TYP(l) -0 THEN GOTO 2575 

2545 READ # 1 7.P . F 1 * ( l . 33 ) 

25S0 IF Fl *(26,33 ) *FNS»( 8,32) THEN 2565 
2*^35 P*P+35 

2556 IF P«075 THEN GOTO 2600 
2560 GOTO 2540 

25H5 WRITE #17.P,F*( 1 ,25>+F2*( 1 .8) .NOENDMARK 
25/0 GOTO 2500 


OHIGlNM- 

of POOR Q uAL11 


>575 WR 
5580 CL 


#1 

#1 


F* ( 1 » 25 ) +F2* ( 1,8) 


mi fl ’- 2 

2600 PRINT M 0 1 6K DIRECTORY FULL " \ M* 1 \CL05E# 1 \ RE 7 URN 

3000 REM COMPUTE AND STORE PROBABILITIES 

3010 PRINT C1*.H1*.B1*.CHR*( 10) ,CHR*( 10) ,CHR*( 10) 

3020 PRINT TAB(25) , "COMPUTE PROBABILITIES" \ 0*1 
3030 FOR R* 1 TO 1000 \ NEXT R 
3040 GOSUB 505 \ REM GET FILE 
3050 IF M«1 THEN GOTO 3120 

3060 PRINT \ PRINT " CAN'T FIND ",F* \ PRINT 

3070 PRINT "Would you like to try again? ", 

3080 INPUT 1 "" ,A* \ IF LEN(A*>«0 THEN GOTO 3080 
3090 IF A * ( 1 » 1 ) * " N " THEN GOTO 3 
3100 IF A * ( 1 » 1 > * " Y " THEN GOTO 40 
3110 GOTO 3080 

3120 GOSUB 990 \ REM LOAD DATA 
3130 PRINT C1*,N1*.B1* 

3140 PRINT "Fliaht name and numbir is " , H 1 * , T* , N 1 * , CHR* <1 0 ) 

3150 INPUT “T/pe in parameter name ",P* \ PRINT 

3160 PRINT "The number of sensor trirrles are ",Y \ PRINT 

3170 INPUT "Tyre in one of the three diffrences ",C7* 

3171 INPUT "Tyre in two of the three diffrences " , D* 

3172 INPUT "Type in thrree of the three diffrences • ,E* 

\ PRINT “DO YOU WANT A HARDCOPY ? " , \ A* * I NCHAR* ( 0 ) 

IF A * * “ Y " THEN S*1 

3193 PRINT C1*.B1* 

3190 IF SOI THEN GOTO 3201 

3191 FILL 65334,84 \ FILL 65532,88 

3192 «#S CHR* ( 27 ) , " COw" , CHR* ( 27 ) , " C2z " 

3200 REM PROBABILITIES 

3201 GOSUB 3210 \ GOTO 3280 

3210 ! #S \ ! #S TAB ( 26 ), “Sensor Redundancy Manaaement" 

3220 ! #S TAB ( 16 ) » "Reror t of Probability Density Function, C ( DF ) " 

2230 1 #S TAB ( 1 2 ) , "Re 1 1 ab 1 1 i t y R, Probalilitx of false Alarm P(AF/DF)" 
3240 1 #S TAB(l7),"For the Sensor Oalue in the Domain of Failure " 

3250 !#S \!#S \'#S TAB( (70-LEN( T* )) /2) , "FLIGHT * ",T* 

3260 '#S "Parameter * ",P*.“ Sensor diffrences are ",C7*," ",D*." and 
3270 '#S \ !#S 
3275 RETURN 
3280 FOR 1*1 TO 
3290 '#S "L-S * 

3300 !#S \ !#S 
3310 !#S TAB (15) 

332u !#S 
3330 P*1 
3340 XI *D( 1 , I ) 

3350 REM COLUMN HEADING FOR TABLE 

3360 '#S TAB ( 10 ) , " PDF " , TAB ( 30 ) , "REL I AB . " . TAB ( 50 ) , " P ( AF/ DF ) " 

3370 '#S T AB ( 10) ,FNS*(53,45) 

3380 REM COMPUTE PROB . 

3380 C*0.3314575#EXP(-(Xl A 2>/2.8974) 

3400 FOR R* 9999 TO 9950 STEP -10 
3410 R 1 * 


31 PRINT 
V 82 S«0 \ 


,D( 1 , I ) , TAB (20) , "L-R 


" » D ( 2 , I ) , TAB ( 40 ) , "S-R * ",D(3.I) 


'P( AF/DF) Table for Sensor Diffrence " . C* . " * " » D ( 1 » I ) 


420 PI 
430 P2 
3440 P3 
3450 P4 
•30 P5 
70 


R/1E4 
C»R 1 
1 - P 1 
1-R1+P1 
1-P3 
f ’30 P5 * 1-R1 

{ 70 A * P1*(P4+P5)/(P1*(P4+P5)+P5*P2) 

3480 '#S T AB ( 10) »C, TABOO) ,R1 .TAB (50) ,A 
3492 P(1,0)*C \ P ( 2 , V ) * A 
3433 0*0+1 
3430 NEXT R 

3500 !#S \ !#S \ P*P+ 1 
351 0 IF P=2 THEN GOTO 3540 
3520 IF P*3 THEN GOTO 3570 


• 53 - 


, E* 


35 JO 
3540 
35 >0 
3580 
35/0 
3580 
30 

- 00 
3605 
3610 
3620 
3630 
3640 
4000 
4010 
4020 
4030 
1035 
4040 
4050 
4055 
4060 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4165 
4 170 
4180 
4190 
4200 
4210 
4220 
4230 
4235 
4 MO 
11 
*..50 
4260 
4270 
4275 
4280 
4290 
4295 
4300 

31*8 

4330 

4340 

4355 
4360 
4365 
4367 
4370 
4380 
4390 
4395 
44 00 
44 10 
44 20 
4430 
4 440 
44 50 
4460 
44 70 
4490 
4195 
4 ’90 
)0 

». .* 1 0 
4520 
4530 
4540 
45l>0 
4555 


' P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE 


" P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE ",E* 


GOTO 3600 
XI* D ( 2 » I ) 

(#S TAB (15) 

GOTO 3360 
XI* D ( 3 » I ) 

•#S TAB (15) 

GOTO 3360 
NEXT I 

IF S*1 THEN A*CALL (65169) 

?B1*,C1*.TAB(20> , ’’CALCULATION COMPLETED" \ PR I NT 
PRINT “ DO YOU WANT A GRAPH ", 

INPUT 1 ” ” , A* \ IF LEN ( At ) *0 THEN GOTO 3630 

IF A*0“Y“ THEN GOTO 90 

REM CREATE GRAPH AND DRAW IT 

PRINT C1*.H1*,B1*,CHR*< 10) ,CHR*< 10) ,CHR*( 10) 
PRINT TAB ( 25 ) , "CREATE GRAPH", Ni* \! 

a * p a i i / c i c q \ 

("SETTING UP PRINTER” ,B1*,CHR*( 10) 

L I NE# 1,132 \FILL 65532,132 \ FILL 65534,(132-4) 


■D» ■ 


! #1 CHR* ( 27 > , ” C 3z ” » CHR* ( 27 ) , ”C4w 
("SORTING DATA " , B 1 * . CHR* ( 10 ) 

REM SORT POINTS IN TERM OF THE X VALUE 
FOR J*t TO 5 

FOR T* J TO 1 5* Y-5 + J STEP 5 
FOR V* T TO 15*Y-5+J STEP 5 
IF P(1,TXP(1,V) THEN GOTO 4140 
T5*P( 1 » T ) \ P< 1 ,T) «P( 1 ,V> \ P ( 1 , V ) «T5 ' 
\ P<2, T) «P(2, V) \ P ( 2 , V ) *T5 > 


\ REM SET PRINTER 


\REM SIZE 


T) 


REM 

REM 


B 1 * > CHR* (10) 
M2 *M 1 \ M4*M3 


THEN 

THEN 

THEN 


Ml 

M2 

M3 

M4 


■P( 1 ,V) 
'P< 1 ,V> 

> P ( 2 , V ) 

> P ( 2 , V ) 


REM 

REM 

REM 

REM 


MIN 

MAX 

MIN 

MAX 


T5*P< 2 
NEXT 
NEXT 
NEXT 

("SETTING GRAPH LIMITS" 

Mi*P( 1 , 1 ) \ M3*P ( 2 , 1 ) \ 

FOR V*1 TO Y* 1 5 
IF Ml >P( 1 ,V> THEN 
IF M2< P ( 1 , V) 

IF M3 >P ( 2 , V ) 

IF M4< P ( 2 , V ) 

NEXT V 

! "CLEARING MEMORY" , B 1 *, CHR* ( 1 0 ) 

REM PLOT IN MEMORY FUNCTION 
V*( 1 , 1 0000 ) *FNS* ( 100,32) 

FOR X*1 TO 10000 STEP 100 \ V* ( X , X+99 ) *S9* 

DEF FNP ( X , Y , C* ) 

V*< (X-l ) * 1 00+ Y ) *C* 

RETUR)' 1 
FNEND 

REM MAKE NUMBER NORAMAL 
("SCALING DATA " , B 1 * , CHR* ( 10 ) 

M5*M2-M 1 \ M6*M4-M3 \ REM TRUE SIZE 

£?1 ,8) = l59( l**?l ,V)-M1 )/M5)*99+l .5) 

P(2»V> *INT ( ( (P(2,V)-M3)/M6>*99+1 .5) 

NEXT 

Po6 L 9»I X VB § RAPH IN mem0RY " »CHR*( 10) 

C*=" "\C*=CHR*(41+J) 

FOR V* J TO ( (Y*15)-5)-5+J STEP 5 
IF P(2,V+5)-P(2,V) *0 THEN H*1 ELSE H*0 
IF H* 1 THEN GOTO 4380 

M * (P( 1 ,V+5)-P( 1 ,V> )/(P(2,V+5)-P(2,V) ) \ REM 
FOR X? * P< ,V) TO P(1 , V+5 ) 

IF H= 1 THEN GOTO 4420 
IF M*0 THEN GOTO 4450 
Y2*P(2,V)-MX2-P( 1 ,V) )/M 
GOTO 4430 
Y2*P ( 2 , V ) 

REM 

Z*FNP(X2, Y2,C* ) 

NEXT 

NEXT 

NEXT 

REM PRINT OUT 

! "PRINTING GRAPH” ,B1* -CHR*( 10) 

S=1 \ REM DEVICE 
GOSUB 3210 \ REM HEADING 
R=MG/10 \ !#S TAB ( 20 ) , \ T*30 
TOR Y7 = M3 TO M4 STEP R 
i * S Y7 , T AB < T ) , 

T*T+10 
NEXT 
! #S\ ! #S 


XCHANG 

XCHANG 


\ NEXT 


SLOPE 




" ,D(2, I ) 

" ,D(3, I ) 

ORIGINAL PAC. 

OF POOR QUALITY 
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*318 ?P r l$oI J0Eiti o 88?o s J§§o 100 

4600 '#S T, TAB (20) , " I " » 

tiiS goto R 4G4o original PAuE rj 

4,630 '#S TAB (20), “I", OF POOR QUALITY 

r »o '#3 o*(x*x+s9> VUMU,r 

<-50 L3«L3+ 1 \ IF L3»ll THEN L3»l 
4660 NEXT 

4665 '#S T+R \!#S\'#S 

4670 !#S \!#S TAB ( 20 >, "CURVE KEY" \ R«9950 
4673 «#S 

4680 FOR X * 5 TO 1 STEP -1 

4690 <#S TAB < 25 ) » R/ 1E4 » " * " , CHR* ( 4 1 ♦X ) , " ",\R«R+10 

4700 NEXT 

4701 <#S 

4710 A«CALL(65169) 

4720 FILL 65532 , 88 \ FILL 65534, (P8-4) 

4725 ! # 1 CHR* ( 27 ) , " C2z" , CHR* ( 27 ) » " COz" 

4710 GOTO 90 

80 vO REM DISPLAY DIRECTORY 
GO JO! B1*,HI*.C1* 

3020 ! “DO YOU WANT A HARDCOPY ", 

6030 A* * I NCHAR* ( 0 ) 

6 ( '40 IF A $ 3 " Y " THEN S»1 ELSE S»0 \ PRINT 
6050 FILL 65532,88 \ FILL 65534 88-4 
6060 OPEN #17.5 , "DIR, "♦STR*(D> 

6070 !#S " FLIGHT AND NUMBER TARGET FILE"\'#S 

6080 IF TYP ( 1 ) «0 THEN GOTO 6120 
6090 READ# 1 , Fl*< 1 ,33) 

6100 PRINT TAB (10), FI $(1,33) 

6110 GOTO 6080 
6120 CLOSE# 1 

6125 A s CALL (65169) \ REM FF 

6130 PRINT "PRESS ANY KEY TO GOTO MEMU ", 

6140 A* = I NCHAR$ ( 0 ) \ GOTO 90 

8090 PRINT Cl *, TAB ( 15) , "PROGRAM TERMINATED" 

8095 PRINT \ PRINT \ PRINT 

mm m STRINGS FUNCTION 

r M0 DEF FNS$ ( L , C ) \ REM L » LENGHT C * ASC CODE MAX LENGHT (255) 
20 S9$* " “ 

1838 l 8 U 9 .lcS 8 .Vc. 


9050 

1898 

9100 
9105 
61 iC 
9 1 20 
9130 
91 33 
9140 
9150 
9160 
9170 
3180 
9100 
9? 00 
3500 
9505 

ISIS 


NEXT 

KM* S9$ 

REM FILE TEST FUNCTION 
REM 0 = NOT FOUND 1,2, 3, 4 IF FOUND 
r NF ( Z9$ ) 

X= 1 

ERRSET 8150, L,E 
X9$ = " , " +STR* ( X ) 

IF FILE ( Z9*+X9* ) <> -1 THEN GOTO 9180 
X* X+l 

IF X< >5 THEN GOTO 9130 

X»0 

ERRSET 
RETURN X 
FNEND 

REM CHECK FOR DIR FILE SET DRIVE DEFAULT 
D=2 \ REM DEFAULT DRIVE 
A = FNF ( "DIR" ) 

IF A = 0 THEN GOTO 9530 


DRIVE LOCATED 


9520 D = A 

9525 GOTO 9550 

9530 CREATE "DIR, 2", 4, 5 

9535 OPEN #1 7.5, “DIR, 2" 

9540 WRITE #1 , & 1 ,NOENDMARK 
9545 CLOSE# 1 
9550 RETURN 

10000 REM TERMINAL CODES 
10010 ei$ = CHR$ ( 7 ) 

10015 Nl$ = CHR*< 126>+CHR*(25) 

) 0020 Cl* = CHR$( 1 26 > +CHR* ( 28 > 
025 HI* = CHR*< 1 2G)+CHR*(31 > 

\ 9 030 Ul* - CHR*( 12S)+CHR*< 12) 
10050 RETURN 
10100 REM DRAW GRID 
10110 FOR X s 1 TO 3 
10120 FOR Y* 1 TO 20 
10130 PRINT FNP* ( ( X- 1 ) *25 , Y+3 > , 
10140 X$*STR*(X) \ Y*=STR* ( Y ) 
10150 PRINT " ( " ,X*(2) , " , " , Y*(2) 


\ REM RING BELL 
\ REM LOW INTENSITY 
\ REM CLEAR SCREEN 
\ REM HIGH INTENSITY 
\ REM UP CURSOR 


l< 


) 


-55 
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10160 NEXT 
10170 NEXT 
10180 RETURN 

10500 REM INPUT MODE PROGRAM 
10510 X5 * 7 \ Y5 *4 \ REM TRUE LOCATION 

10520 Y6*0 \ YG«0 \ REM END LOCATION 

530 A*=FNS*< 11.32) 

4 540 PRINT FNP*(X5.YS). 

10550 IF Y6*57 AND Y6*23 THEN GOTO 10690 
10560 GOSUe 11410 
1 0570 IF Y»0 THEN GOTO 10G90 
10580 05* a I NCKAR* ( 0 ) \ 05=ASC(Q5*> 

10600 IF 05 >47 AND 05<58 OR 05*46 OR C5*69 OR 05*45 THEN GOTO 10G20 
10610 GOTO 10720 
10G20 N5* = " " 

10630 GOSUe 11140 \ REM INPUT # 

10G40 X6*X5 \ YG*Y5 

IOG50 X5*X5+25 

106G0 IF X5<>82 THEN 10540 

10670 X5*7 \ Y5*Y5+ 1 

10675 IF Y5<> 24 THEN 10540 

10678 X5-57 \Y5*23 \ GOTO 10540 

10G90 05** INCHAR* ( 0 ) \ Q5*ASC(Q5*) 

10700 IF 05=67 OR 05*68 OR 05*83 OR 05*76 OR 05*82 OR 05*85 THEN GOTO 10720 

10710 GOTO 10690 

10720 IF Q5*v >"S" THEN 10740 

10730 IF X6 *57 AND YG>3 THEN RETURN 

10740 IF 05*'. > M U" THEN GOTO 10760 

10750 IF Y5<>4 THEN Y5*Y5-1 \ GOTO 10540 

10760 IF 05* < S >"D" THEN 10820 

10770 R*X5 \ C»Y5 \ REM TEMP 

10780 Y5= Y5+ 1 

10781 X-X6+25 \ Y*Y6 \ IF X< >82 THEN 10783 

10782 X*7 \ Y* Y6 + 1 

10783 IF Y< >24 THEN GOTO 10784 \ X*X6 \ Y*YG 

10784 IF X »X5 AND Y*Y5 THEN GOTO 10540 
10790 GOSUB 11420 

10800 IF Y*0 THEN GOTO 10540 
10810 X5=R \ V5«C 
10815 GOTO 10540 
1( ’»820 IF Q5*< >"L 1 THEN 10850 
330 IF X5< >7 THEN X5 = X5-25 
J840 GOTO 10540 
10850 IF Q5*<>"R M THEN GOTO 10930 
10860 IF X 5 = 5 7 THEN GOTO 10540 
10865 GOSUB 11420 \ IF Y*1 THEN GOTO 10540 
10870 X5*X5+25 

10880 IF Y50Y6 THEN GOTO 10910 
10890 IF X5<*X6+25 THEN GOTO 10540 
10900 X5*X5-25 

rnt ?P T 85i??«- THEN GOTO 10540 

10935 GOSUB 11420 \ IF Y*1 THEN GOTO 1054t. 

10940 PRINT FNP*(X5. Y5> .A*. 

1 0950 PRINT FNP*(X5,Y5). 

10960 GOSUB 11110 
10970 GOTO 10540 
11110 REM INPUT NUMBER 

mss m * I NCHAR* ( 0 ) \G5 = ASC ( 05* ) 

REM GET 

IF 05=13 THEN GOTO 11280 
IF 05=8 THEN GOTO 11220 
IF LEN < N5* ) = 1 0 THEN GOTO 11130 
N5*=N5*+G5* 

PRINT 05*. 

GOTO 11130 

IF N5* = " " THEN GOTO 11130 
IF LEN ( N5* ) < > l THEN GOTO 11250 
N5* = " " \ GOTO 11260 
N5* =N5* < 1 . LEN ( N5* ) - 1 ) 

PRINT 05*. 

GOTO 11130 
ERRSET 11310. L.E 
D ( ( X5-7 ) /25+ 1 , Y5-3) =0AL(N5* ) 

ERRSET \ RETURN 

PRINT FNP* ( X5 , Y5 ) . " ERROR ".Bl*. 

FOR Z=1 TO 500 \ NEXT Z 
PRINT FNP* ( X5 . Y5 ) . A* . FNP* < X5 . Y5 ) . 

ERRSET 
GOTO 11120 

REM POSITION LOCAT SIZE 

Y = 1 \ REM Y= 1 IF TRUE >END ELSE 0 -56- 

IF X5*X6 AND Y5=Y6 THEN Y=0 




11140 
11190 
11200 
11201 
1 1 205 
1 120G 
11210 
11220 
11230 
11240 
11250 
1 1260 
11270 
11280 
; t 790 
300 
\ x310 
11320 
11330 
11340 
1 1 350 
11410 
1 1 420 
11425 


V 

11430 IF X5 > «X6 AND Y3<Y6 THEN Y»0 
114 40 IF X3< X6 AND Y5C-Y6 THEN Y«0 
11450 RETURN „ 

30000 REM FNT* FUNCTION TO STRIP OFF TRAILING SPACES 
3«>003 REM T* CONTAINS STRIPPED WORD L IS THE NEW LENGHT 
30005 DEF FNT * ( T* ) 

">07 X»LEN ( T* ) 

3,010 FOR L»X TO 1 STEP -1 

30015 IF T* ( L * L ) < > " M THEN EXIT 30035 

30020 NEXT 

30023 T*»" “ \ L«0 

30030 GOTO 30040 

30035 T*» T* ( 1 , L > 

30040 RETURN T* 

30045 FNEND 

30100 REM PLOT TO SCREEN (TERMINAL) 

30110 FEF F N P * ( X * Y ) 

30120 RETURN CHR* ( 1 26 > +CHR* ( 1 7 > +CHR* ( X ) *CHR* < Y > 

30130 FNEND 

§ 0000 REM OAR I ABLE TABLE 
0005 2 TAB STOPS (TEMP) 

GOO 10 X FOR NEXT LOOPS (TEMP) 

60015 S9* STRING* FUNC . (SET) 

60020 Bl*>Nl*>ClSrHl* TERMINAL CODES (SET) 

60025 29* FILE TEST FNF (TEMP) 

60030 X9* TEMP VARIABLE (TEMP) 

60035 A TEMP VARIABLE (TEMP) 

60040 A* TEMP VARIABLE (TEMP) 

G0045 F* INPUTED FILE NAME 
60050 FI* DISK FILE NAME 

§ 0055 F2* TRUE FILE NAME 
00GO T* STRIPPED WORD 
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■ Sensor Redundancy Management • 

This la a program to allow the input o* eeneor data to the dati 
♦lie ♦or proceeding by the report generating program. 

NASA Grant NAG 4-* 

June 1401 

Prank Scalro. Ph.D. 


— DISPLAY HEADING ON THE SCREEN — 


12 CLSi CLEAR lOOO 'met aside string space 

13 PRINT 968. CHFrS ( 23) t "SENSOR REDUNDANCY MANAGEMENT" 

14 PRINT STRZNGS(32. “•"> 

15 PRINT 9324 9 "Probed i 1 i ty o* False Alarm" 

IS PRINT TAB (S) | "Data Entry Program- 

17 PRINT 9432 # "»y F. Scalzo. Ph.D." 

20 * 

21 ON ERROR GOTO 2000 

22 OPEN M* # l t -SFILE1" 

23 CLOSE 1 

24 PRINT "There is already data on file" 

23 INPUT "DO YOU NANI TO REPLACE I T " | AS 
2A IF LEFTS (AS. 1 ) <> "Y" THEN RUN -MENU- 
27 CLS 

30 INPUT "Enter the number oi Sensor Pairs "|N 
3S DIM S (N. 2) 

40 FOR J-l TO N 

43 PRINT 9143, CHRS (31 > 'Erase to the end o4 screen 

30 PRINT 9143. "Enter the Data For Pair S"|J 

33 PRINT 9323.1 

SO INPUT "READING bl "iS<J,H 

S3 PRINT 9387 # | 

70 INPUT -READING 42 -|3<J,2> 

•O NEXT J 
40 PRINT 

100 INPUT -Would you like to see the data "|AS 
110 IF LEFTS < AS . 1 > <> "V THEN GOTO 140 
113 CLS 

120 FOR J-l TO N 

123 IF INT(J/»2> - J/ 12 THEN GO SUB lOOO 

130 PRINT "Pair "| J, -Reading #1 - 1 3 < J , 1 ) , -Readi ng S2"|3<J,1) 

140 NEXT J 
130 CLS 

160 INPUT -Do You Have Any Corrections "tAS 
170 IF LEFTS(AS,1) 0"Y" THEN GOTO 240 
190 INPUT -Which Pair " j J 
140 INPUT "Reenter Reading SI "tS<J,l> 

200 INPUT -Reenter Readirg #2 "lS(J.2> 

210 INPUT -Any Other Corrections ‘iAS 

220 IF LEFTS < AS, 1) - -V" THEN CL3« GOTO 190 

230 CLS 

240 INPUT -Do You Want a Hard Co©y of the Data "|AS 
230 IF LEFTS < AS. 1 ) <> -V" THEN GOTO 300 
260 FOR J-l TO N 

270 LPRINT -Pair-j J, "Reading 41 - i S (J. 1 ), "Reading #2"tS<J.l> 

290 NEXT J 

300 OPEN -0", 1. ’3FILE1" 

310 PRINT SI, N| 

320 FOR J-l TO N 

330 PRINT SI, S(J,1>«S(J,2>| 

340 NEXT J 
330 CLOSE t 
400 RUN -MENU" 
son end 

1UO0 PRINT 9400. " (press any key to continue>"l 
101O IF INKEYS--- THEN GOTO 1010 ELSE CLS 
1020 RETURN 

2000 * ftttt* FILE DOES NOT EXIST — OK TO CONTINUE tSSift 
2010 RE3UME 27 
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.999 
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.0*99112 

.9*43 
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.0*9*312 
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.0*9*312 

.9*42 
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.9*41 
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.0*991X2 

.9*4 
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.0*9*112 

.99*9 
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.911702 

.0*991X2 

.9919 

.9141* 

.0*991X2 
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.0*991X2 

. 9*93 

.*12X42 
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.9914 

.911X21 

.0*9*112 

.••12 
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.0*9*112 

. 9*32 

.9490** 

.0*991X2 

.9*11 

.94*0** 

.0*993X2 

.991 

.9490* 


0fM0/O0t 

090 
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.032X441 
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.99*4 
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.0121441 
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.932342 

.032X441 
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. 949722 
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.***2 

944914 

.032X441 

.99*1 

.944123 

.032X441 

.99* 

.94X333 

.032X441 

. 9*7* 

. 97*394 

.012X441 

.9*7* 

.912*31 

.032X441 

. 9*77 

.973173 

.0321441 

• 9*7* 

.930413 

.032X441 

.9*73 

.*27714 

.032X441 

.9*74 

.*=3034 

.03214*1 

.9*73 

.. *2234* 

.0321441 

.9*72 

.91*714 

.032X441 

.**71 

.91707* 

.032X441 

.**7 

.914439 

.0321*41 

.***9 

.911*4* 

.032X441 

.9*4* 

.90*234 

.0321441 

.9*47 

.904474 

.0321441 

.9*44 

.904X0* 

.0321441 

.**43 

.90X537 

.0321441 

.9*44 

.***019 

.0321441 

.9*43 

. *944*3 

.032X441 

.9*42 

.*939*3 

.0321441 

.9*41 

.*9X4** 

.0321441 

.•*4 

.**9004 

.0321441 

.9*-* 

. **4333 

.0321441 

.9*30 

.*•4077 

.0321*41 

• 9*37 

.*•1433 

.03214*1 

.**34 

.*79202 

.03214*1 

.9*33 
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.*7437* 
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.**33 

.*719*4 
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.9*32 

.*4*404 
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.9*31 

.*4724 
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.*44**4 
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PROGRAM SENTTAT 
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FTN 4. 2*75060 


( 


{ 


PROGRAM SENSTAT (INPUT, OUTPUT, TAPE5*INPUT,TAPE6*OUTPUT, 

♦ TAPE1,TAPE2,TAPE3,TAPE4) 

C 

COMMON/STATOAT/STATMIN,STATMAX,MEAN, VAR,OEV,MEANSO,RMS , 

* ♦ SKEW, KURT, XLOW,XHIGH,NSGOOO,NSTOTAL 

COMMON/MISTOAT/HISTMIN,HISTMAX *nsegs,tally,nhgooo,nhtotal 
COMM 3N/PLOTOAT/XSCALE (4) , YSCAlE (4) , H ST AR T, HSTOP, TI TLE ( 6) , IOPLO T 
REAL MEAN, MEANSO, KURT 
INTEGER T ALLY (10 0) 

01 ME >< SI ON OATA (100 ) , ITIME (6) , MESSAGE ( 3) 

C 

REA0(5, 101 TITLE 
10 FORMA T ( 5 A 1 0 ) 

IF (EOF (5) .NE.O)STOP -NO TITLE CARO" 

C 

REAO( 5, 10 ) HST ART , HSTOP 

IF (EOF (51 .NE.OISTOP “NO TIME CARO" 

BACKSPACE 5 
C 

RE AO( 5,20) ITIME 
20 FORMAT (2(3 (12, IX) , IX) ) 

ST ART sFLOAT ( 10 0 0 * ( 3600*ITIME(1) *60*ITIME(2) ♦ IT IME ( 3 ) ) ) 

STOP* FLOAT ( 10 0 0* ( 3 600* IT I ME (4) *60* IT I ME (5) ♦ IT I ME (6) )) 

C 

RE*Q|5|3Q> IOPLOT,lNOEXA, I NOEXB , NPARA MS 
30 FORMAT | A10. 3110) 

IF |EQF ( 5 1 » NE.OISTOP "NO JOB CAROS" 

C ALL PLOTS (0,0 ,4) 

CALL FACTORI. 7671) 

CALL PLOT I 0 « . 1 1 ♦ »- 3) 

GO TO 50 
C 

40 RE AO( 5,30) IOPLOT,INOEXA,INOEXB,NPARAMS 
IF (EOF ( 5 ) • NE, 0 ) GO TO <>93 
C 

50 IF (INOEXA.GT.O .ANO.INOEXA .LE.NPARAMS. ANO. 

* IHOEXB.GE.C.ANO.INOEXB.lE.NPARAMS)GO TO 55 
CALL REMARK(-ILLEGAL JOB CARO") 

• GO TO 40 
55 REWINO 1 

RE wi NO 2 ORIGINAL HALE IS 

N*0 Qp poOR QUALITY 

MESSAGE(l) * "WRITING 2" 

MESSAGE (2) =IOPLOT 

MESSAGE ( 3) =0 

CALL RFMAPK (MESSAGE) 

C 

60 REAO( 1 ) Tl,T2,TMILLIt (OAT A (I) *1*1, NPARA MS) 
IF(EOF(1).NE.O.OR.THILLI.GT.STOP) GO TO 70 
IF (TMILLI.LT, START) GO TO 60 
A * OAT A ( INOEXA) 

8 = 0 . < 

IF (IN0EX9.NE.C) B = D AT A (INOEXB) 

C= A-3 
WR IT£ ( 2) C 
N=N*1 
GO TO 60 
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PRO GR AM SENSTAT 


73/74 0® T* 1 


FTN 4# 2 ♦ 75 06 


C 

70 EMOFILE 2 

IF(N.EO.O)STOP -no OATA within time slice- 
ST ATM I N*-l E99 
ST ATI A X* ♦ 1 E99 
* MESSAGE ( 1) ***S T ATIST Z~ 

CALL REMARK (MESSAGE) 

CALL STATIST (2) 

M I STM I N*XL OW .. 

HISTMAX*XMTGH ORIGINAL * kU-14 

MSECS *100 OF POOR QUALITY 

CALL HSTGR AM ( 2 ) 

IF(IN0FX9.EQ.C )GO TO 40 
C 

REHINO 2 
REWINO 3 

HISTMInsHISTMA X*0. 

MESSAGE CD* -WRITING 3" 

CALL REMARK (MESSAGE) 

C 

80 RE AOl 2 ) C 

JF (EOF ( 2) • NE. 0 ) GO TO 90 
X*ABS CC-MEAN) 

HIST1AX = AMAX1 ( hi st MA X « X) 

HR ITS I 31 X 
GQ H SO 
C 

90 ENOFJLE 3 
NSEGS* 100 
CALL HSTGR AM( 3 ) 

GO TD 40 
C 

990 CALL PLOT(20.*0. ,999) 

C 

STOP -ENO OF JOB" 



UBROUTINE STATIST 73/74 0>T«l , , rl n. cerT 13 FTN 4, 2^75060 08 

origin* 1 - ITV 

OF POOR QUALITY 

SUBROUTINE STATIST (FILE* 

COMMON /ST ATOAT/XMlN, XMAX, MEAN, VAR, OEV, ME ANSO, RMS, SKEW, KURT , 

♦ XLOH,XHICH,NGOOO,NTOTAL 
COMMON/PLOTOAT/XSCALE (4» , Y SC ALE (4 I , H ST AR T, HSTOP, TITLE ( 6) .IOPLOT 

• REAL MEA N, MEAN SQ, KURT, N 

INTESEP FILE 
C 

ME an* VAR* OEV* MEANS Q* RMS* SKEW*KUR T*SUMX*SUMX2*SUMX3*SUM X4*0 . 

XLOM* ♦ 1E99 
XH IG^ *- IE 99 

ngooo*ntotal*o 
rewino file 

C 

10 REAO(FILS)X 

IF (EOF (FILE) • NE, 0) GO TO 20 
NTOTAL*NTOTAL*l 

IF(X.LT.XMIN.0R.X.GT,XMAX)G0 T9 10 
NGOOO * NGOOOH 
SUMX* SUMXf X 
X 2*X* X 

SUMX2*SUMX2*X2 
X 3 *X* X 2 

SUMX3*SUMX3*X3 
X4 = X*X3 

SUMX4*SUMX4*X4 
XL0« s A MINI ( XL O W , X ) 

XHJGi* AMAX1 (XH IGH, X) 

GO TO 10 
C 

20 IF (NIOOO.EQ.O) STOP -SUBROUTINE STATIST HAS NO OATA~ 

NsFLO A T ( NGOOO* 

ME AN*SUMX/N 
VAP=SUMX2/N-ME AN** 2 

IF(VAR.LE.C.»ST0P -SUBROUTINE STATIST SHOWS ZERO VAR I A NCE” 
OEV=SQPT(VAR) 

ME AN30=3UMX2/N 
RMS=SQPT (MEANSQ) 

SKEW* (SUMX3/N-3, *MEAN*ME ANS0*2 ,*MEAN**3) / (OEV* VAR) 

KURT* I SUMX4/N-4. *MEA N* SU MX3/ N* 6 . * ME A N* *2 *SUMX2/ N-3 . *ME AN**4) / 

♦ ( VAR* VAR) 

WR ITE ( f>, 3C ) TITLE, FILE, IO PL OT,N TOTAL, XM IN, XMAX, NGOOO, XL OW, XHIGH, 

♦ MEAN, VAR, DEV, ME ANSQ , RMS, S KEW, KURT , HST ART, HS TOP 
30 FORMAT (1H1////- RESULTS OF SUBROUTINE STATIST**//// 

1 - TITLE”. T28,6A1C//" OATA FROM FILE NO. - ♦ T 28 , 1 2/ / 

2 - PLOT I 0-.T28, A 10//- TOTAL NO. OF SAMPLES”, T2 8, 16// 

3 ” ALLOWABLE MI NI MUM- , T2 8, G12 . 6// 

4 - ALLOWABLE MAXIMUM-, T28.G12. 6// 

5 - SAMPLES ALL0WE0",T28,I6// 

6 ” OBSERVEO M IN I MUM" , T2 8, G 1 2. 6/ / 

7 ” OBSERVEO MAXIMUM", T28, G12. 6// 

8 " MEAN VALUE",T28, G12.6//" V AR I ANCE", T 28 , G 12. 6// 

| 9 ” STANOARO DEVIATION", T28,G12. 6//" MEAN SQU ARE O” , T 28 , Gl 2 . 6// 

A - ROOT MEAN SQUARE 0" , T 28 , G 1 2 . 6 //" SKEWNESS”, T28, G12. 6// 

9 ” KURT0SIS",T28,G12.6//" START T IME", T28 , A 10// 

C ” STOP TIME",T28,A1Q) 

RETURN 

ENO 


- 63 - 


i> 


I0ROUTINE HSTGRAM 


7 3 / 7 <♦ OPT * 1 


OF POOH QUALITY 


FTN 4.2»75060 


SUBROUT I NE HSTGRAM (FILE) 

C 

COMMON/HI STOAT /XMIN.XMAX, NSEGS »?ALLT»NG000,NT0TAL 
COMMON/PLOTOAT/XSCALE (4) ,YSCALE<4) ,HSTART,HSTOP»TITLE(6) , I0PL0T 
CO MMON/OAYTIME/OAY .HRMlNSC 

* INTEGER FILE,TALLY(100)»TOTAL<100» ,REHAINC100I .PASSNO 
01 ►CNS ION FROM C 100 ). TO (100). MESSAGE C 4) 

C 

MESSAGE ( 4) *PASSNO* 0 
PERCENT *0. 95 
NCONSE C* 3 
C 

1 PASSNO=PASSNO* t 

ENCODE (2 7, 5. MESS AGE) FILE,PASSNO, IOPLOT 
5 FORMAT ("HSTGRAM Ml," RUN M ,I2,1X,A10) 

CALL REM APK (ME SS AGE) 

C 

00 10 1*1, NSEGS 
10 T ALLY ( I) *0 
C 

REMINO FILE 
ngooi*ntotal*o 

' OEIU = 1.00 000 1 MXMAX-XMIN) /FLOAT (NSEGS) 

c 

20 RE AQ( FILE) X 

IF (E)F (FIVE) ,NE. 0) GO TO 40 
NTOTAL*NTOTAL* 1 
IF (X, L T » XHlN) GO TO 20 
TEST* X M JN 
C 

00 30 1=1, NSEGS 
TEST*TEST»OELTA 
IF (v. z i . TEST) GO TO 30 
TALL* ( I) =TALLY (I) M 
NG009*NG000U 
GO TO 20 
30 CONTINUE 
C 

GO TO 20 
C 

40 CALL DATE(OAY) 

CALL T IME (HRMI NSC) 

WRIT; » 6,50) TITLE .FILE, IOPLOT,NTOTAL, XMIN,X MAX, NGOOO, NSEGS, 

♦ OAY.HRMINSC 

50 FOPHA T ( 1H1//" RESULTS CF SUBROUTINE HSTGRAM**/// 

1 ** TITLE". T28.6A10//** OATA FROH FILE NO. ** , T 28 , 1 2/ / 

2 - PLOT I 0",T28, A10//" TOTAL NO. OF SAMPLES", T2 8, 16// 

3 ** ALLOWABLE MI NI MUM" , T2 8 , G 1 2 . 6 // 

4 " ALLOWABLE MA X I MUM" , T 28 , G 1 2 . 6 // 

5 " SAMPLES ALL0WE0",T28,I6// 

6 " NO. OF SEGMENTS", T28, 13// 

7 ** OATE/T IME",T28,A10 ,A9//> 

C 

TOTAL ( 1) sTALLY (1) 

REMAIN (1) *NGOOO- TALLY Cl) 

FROMC 1 ) s X = X M I N 



U0ROUTINE HSTGRAM 
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0 


C 

c 


c 


c 

c 


c 

c 

c 

c 


c 


c 

c 


00 60 1*2, MSECS 

TOTAL ( I) ’TOTAL C I-t ) »T ALL Y C I) 

REMAIN Hi «NGOO 0- TOTAL C I) 

.60 FROMC I ) *T0 ( 1-1 ) *X*X’OELT A 


ORIGINAL PAGE 13 
OF POOR QUALITY 


’ TO (NSEGS) «XMAX 

WRIT- 16,70 (I, PROM Cl I , TO flit TALLY CD » TOTAL III , REMAIN (I ) , I* t, NSEGS) 
70 FOPMATC- SEC. NO FROM" , T29, "TO", T 3 8, **S AMPL ES TOTAL**, 

♦ T5 7,**RE MA INING"// 100 ( 2H ,I3,4X,2(G12.6«2X) ,3(16,4X1/) ) 


call plhist 

OLOMI N’XNIN 
OL OMA X * X M A X 

MlNR-Q* I NT (PERCENT ’FLOAT (NCOOO) ) 

IF(FILE»EO. 3) CO TO 110 
IFCPASSNO.EO. DCALL PLNORM 
I M I N* NCONSEC’i 
Ml NO* 0= (NGOOO’MINRFC) /2 

00 100 I*IMIN, MSECS 

IFCREMAINC I-l) .LT. MINREQ)G0 TO 110 

DO 03 J* 1 , NCONSEC 

K= J*4 

|F (TALLY CXI .NE 'CO TO 100 
90 CONTINUE 

XMIN*FPOM< I) 

100 CONTINUE 

110 I s NSEGS- NCONSEC 

120 IF (TOT AL C I ) .LT.M INREO) GO TO 150 

00 130 J*l, NCONSEC 
<= I*J 

IFCTALLY(K) .NE. 0)CO TO 140 
130 CONTINUE 

XMAXsTOCI) 

140 1=1-1 

IFCI.GF.HGO TO 120 

150 IF(XHlN.NE.OtOMIN.OR.XMAX.NE.OLOMAX) GO TO 1 
RETURN 
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SUBROUTINE PIHIST 
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SUBROUTINE PLHIST 


original page W 

OF POOR QUALITY 


FTN 4.2* 75060 


< 


C 

C 

C 

c 


c 

c 


10 


COMMON /H IS TOAT /HIS TMlN,H I STM AX, NSEGS, TALLY, NHGOOO.NM TOTAL 
C0HH3 N/pl 0T0AT/XSCALE<4» , VSCALET4I ,HSTART,MSTOP, TI TL E t 6 1 , I OPLO T 
COMMON/O AY TIME/O AY, HRMINSC 
INTEGER T ALLY (13 0) 

01 HEN S I ON T0PLINEI4) 

OATA TCPLINE/~MINIMUM MAXIMUM VALIO 08", "S"/ 

CALL PLOT (11. 43,0, ,-3» 

XSCALE (1 l*HISTMLN 
XSCALE ( 2 I *HIST MAX 
call SCALE(XSCALE,10.,2,H 

CALL AXIS(0.,0 .,“SCAL AR VALUE"»-12 ,1 0. ,270 . , XSCALE ( 3 1 , XSCAl E ( 4 )) 
YM AX* 0 . 

OO 10 1*1, NSEGS 
TESTsFLOAT (TALLY (I ) ) 

YMAYsAMAXl (YMAX.TEST) 

yscale ( i) *c. 

YSCAl e ( ?) * ymax 

call STALE (YSCALE, A. ,2,1 ) 

CALL AY90(0., G ., "SAMPLES PE® INTERVAL' 


,20, A. , 0 . , 0. .YSCALE (4) ) 


X* (Hf S TMlN-XSCALE ( 3) ) /XSCALE ( 4) 
CALL ®LOT(0.,-1.*X,3) 


OELTAXsABS(MlSTMAX-HISTMINI/FLOAT(NSEGS> 

OO 20 1*1, NSEGS 

YsFLOAT (TALLY ( I) )/ YSCALE ( 4) 

CALL PLOT ( Y,-l ,*X, 2) 

X*x*0ELTAX/XSCALE( 4) 

20 CALL PL0T(Y,-1.*X,2) 

CALL ®LOT( 0,-1 ,*X, 2) 

CALL SYMBOL < 0 . , l , 5 , . 1 5 , T I TLE , 0 . , 60 ) 

CALL SYMBOLIC. ,-5.87,. 12, TOPLINE, 2 70., 31) 
CALL NUMBER (8. 7 5, -5. 97,. 12, HIS THIN, 2 7 0 ., 3) 
CALL NUMBER (8. 75,-7.35,. 12, HIS TH AX, 270. ,3) 
HGOOO = FLOAT(NHGOOO) 

CALL NUM8ER(8.75,-8.50,.12,HGOOO ,270.,-l) 

CALL SYMBOL (8. 7 5, -9. 3 5,. 1 2 ,"/", 270 . , 1 ) 
HTOTAL*FLOAT (NHTOTAL) 

CALL NUMBER(8. 75,-9.50,. 12, HTOTAL ,270. ,-l) 
CALL SYMBOL (8. 0,-.5,.12, IOPLOT,270., 10) 

CALL SYMBOL(8. 0,-2.,. 12, HSTART ,270., 8) 

CALL SYMBOL (8. 0 ♦ -3 . 12 , . 1 2 , "TO" , 270 . , 2 ) 

CALL SYMBOL (8. 0, -3 .50 , .12, HSTOP, 270. , 8) 

CALL SYMBOL (8. 0,-7.75,.t2,OAY,270. ,10) 

CALL SYMBOL (8. 0 , -9 . 0 0 , . 1 2 , HR MI NSC, 27 0 . , 9) 

RETURN 

END 
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[UBROUTINE PLNORM 73/74 OPT«t FTN 4.2*75060 

; ORIGINAL J* 

tr OF POOH QUALITY 

SUBROUTINE PLNORM 

c 

CO MHO N/ ST A T OA T/ STA TM IN, ST ATMA X , ME AN , V AR ,OE V, ME ANSO, R MS . 

* ♦ SKEW, KURT .XLON.XHICH.NSGOOO.NSTOTAL 

COM 13 N/HISTOAT/HISThIN.MISTmAX.NSECS, TALLY, NHC 00 O.NH TOTAL 
COMMON /PLOT OAT /XSCALE(4) ,Y SCALE (41 , HST AR T, MSTOP. TITLE I 6) ,I0PL0T 

real mean,meansq,kurt 

INTESEP TALLYdOul 

01 MEN SI ON ABSCISS ( 302) ,0 PON ATE ( 302 I, TOPLINE (6) 

C 

OATA TOPLINE/-HEA7 SOUAR","E RMS VARIANCE". “ MEAN VA 

♦ "LLE STO"," 0EW 

C 

x*xs:ale(3) 

XINC*APS(XSCALE(4) /2 9. 9) 

COE F F * FLOAT (NS GOOOJ* (XHIGH-XLOW) /FLOAT (NSECS)/ (2. 50 66* OEV) 

C 

00 10 1*1,300 
A9SCI SS ( I) *X 

EXPOV*-. 5* ( ( (X-MEAN) /OEV) *’2) 

ORONATf (I)sO, 

IF(E<P0N.GE.-6 75 ..AN0.EXP0N.LE.74l.)0R0NATE(I) *COEFF *E XP ( E X PON ) 
10 X* X*X I NC 
C 

ABSCISE(301)sXSCALE(3) 

ORONA TE ( 301) *YSCALE(3) 

ABSCl S? ( 302) *-XSCALE (4) 

OPDNATE I 302) *YSCALE(4) 

C 

CALL L l NE (ORONATE, ABSC ISS , 30 0 , 1 , 0 , 0 ) 

C 

CALL SYMBOL ( 9 . , l . , . I 2 , TOPL INE , 27 0 . , 5 4 ) 

CALL NUM9ER(8- 75, *0.60 ,. 12.MEANSQ ,270 . ,3) 

CALL NUMBER(8. 75.-0.70 ,. 12, RMS ,270. ,3) 

CALL NUMBER(8. 75,-1. 85,« 12, VAR ,270., 3) 

CALL NUM8ER(8.75,-3.25,. 12, MEAN ,270. ,3) 

CALL NUMBER(8. 75,-4. 75,. 12.0EV ,270. ,3) 

C 

RETURN 



PROGRAM NOR HP I T 


71/74 0PT*1 


' V 






( 



ORIGINAL FACC 13 FTN4.2»75 060 

Of POOR QUALITY 

p°nc?4 NOPHfiT ( input ♦output .tapes * input ,tapeb*output, 

♦ TAPE1,TA0E2,TAPE3,TAPE4» 

C 

C0HM0M/PTAT0AT/STATMIN,STATHAX.K£AN« VAP.OEV* HE ANSOtRMS . 

♦ S<E W, KURT, XLOW,X HIGH, NSGOOO, NSTOT AL 
C0HM3N/msT0AT/HlSTHIN,HISTH4X,N$EGS, T ALLX, NHGOOO, NHTOTAl 
COHM-JN/OLOTOAT/X^CALE (4) « YSC ALE (4) , HS T AR T , HS TO®, T I TL c ( 61 , IOPLOT 
REAL Hf 4 n,m£AnS3,KURT 

INTESP“ TAUMAOO) 

01 hr I ON DATA (130), I TIME <6) , HESS AGE 131 
C 

REAO( K , l 0 ) TITLE 
10 P0 c hA T M A 1 C ) 

I F ( R} R ( 5 ) , N£ , C ) STOP "NO TITLE CARO” 

C 

R£ 7-D( 5 » l P I «START ,HSTOP 

IPt r 3*'l5) .NE.OSTOP "NO TIHP CARO" 

PAC;5 e fCP 5 

C 

°F AD( C , ?C ) I T I HE 
20 POPHftT C?(3fI2, 1X1 , IX) ) 

S T APT • FLOAT ( 10 00 *<3600*1 TIME 11) ^60* ITI HE C2) ♦ITIHEC3H) 

S T O®* FLOAT ( 13C0*(J6Q0*ITIMF(4) *60*ITIHE (5 ) ♦ITIhE (6 ) 1 ) 

c 

REAO( 5 ♦ JO) IO°LOT .INOEXA, INO£XP,NPARAHS,TOL 
30 p OPHH |AU , 3110, F20.6) 

JFIEOFm .ME.CISTOP "NO JOB CAROS*' 

CALL PLOT? IC.G .4) 

CALL FACTO® (. 7 871 ) 

CALL Pl^T(C..ll. .-3) 

GO n fc j 

C 

4 0 ®E AO I *5. 30 I IO°LOT,lNnEXA, I NOE Xfl , NPAR A MS , TOL 

IF (COF (S) .NE.C ) GO TO 990 
C 

5 0 IFflNOf XA.GT.C . ANO.INOEXA.LE.NPARAHS, A NO, 

♦ INOF XT. CT, 0 . ANO . IHOEXB.LE . NPARAHS ) GO TO 55 
CALL P‘ MARK ("ILLEGAL JOB CARO**) 

GO T) 4J 
55 OPWINO 1 
rewind ? 

H- C 

ME*FA GF (1) * ,,, .IP ITING 2** 

HE ?SA G r (?) sIOPLOT 

HEFSAGf ( J) s0 

CALL Pf HtFK C1ESSAGF) 

C 

6 0 OFAO( 1 ) T1,T2.THILL I, (OAT A (I) , I~ 1 , NPAR A HS) 

IF (Fri r ( 1 ) . NF ,C .OR. THILLI . GT. STOP) GO TO 70 
IF (THILLI.LT, START) GO TO 60 
C = DAT A ( INDFXA) -0A T A ( I NOE Xfl ) 

W»IT: (2)0 
N*NU 

go n f 0 

C 

7 0 ENOFiLf 2 

IF (N. EC , o STOP "NO DATA WITHIN TIME SLICE" 
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PPOGOAM NOP HP IT 


73/74 


0PT*i 


PTN 4 . ?♦ 75060 


ST£TNTH«-1E99 
STATnAX*» 1E99 
HE SSA GE < i I * "ST atist 
CALL °f '1APP ('if 9SAGE) 

CALL ST ATIST (? I 
* XLOW*hI5THINs-i,?*TOL 
XM IG-o HIST MAX* M .2 *TOL 
SSPG5* 10? 

CALL h$TGPAm<2) 
NSrOOH*NMf,OOT 
CALL PLNOPM 
C 

pewmo ? 

XLCMsM JTT HI N = HP A ‘1-3.* nev 
XHTGHsHiSTHAXxHE AN*3.*OE V 
CALL T r .PAM<2 ) 
NSc.ODO^NHr.oc*) 

CALL 

go ro 40 

C 

99C CALL P| 01 (?0. • ). .999) 

C 

STOP ••‘■NO OP JO?** 


original rr.cz rr 

OF POOR QUALflV 



SUBROUTINE STATIST 7 3/74 OPT«l ORIGINAL PA«E » FTN 4.2*75060 

I . , OF POOR QUALITY 


r 

< 


f 



c 


c 


c 


SU^DUTJN? STATIST (FILE) 

COMMON/ STATOAT/XMlN, XMAX, MEAN, VAR, OEV. ME ANSO. e MS, SKEW, KU»T, 

♦ XI.0W, XHICM, NGOOO* NTOT AL 

COMMI^/^rTTAT /XSCAtE (4) .YSCALE 14) .HSTART,HSTOP. TITlF(6l » I OPLOT 
PEAL MF AN,M;ANSO,KURT,N 
INTEOEF FILE 

MEAN*VAR*nEVsMEANSO«RMS*SKEW *KUR T» SU MX * SUM X 2* SUM X 3 » SI 1“ X4 =0 . 

XLOW* ♦ IE RQ 

*Hlf,Hs-lf9R 

NGCOO*NTOT AL«C 

OF VII SO r I Lt 

19 REATIFIL ; >X 

IF (F3r ( FILE) . NE. 0) GO TO 20 
NTOnL s NTOTAl» i 

IFCX.L t .xmIN.C’.X.GT.XMAX)GO TO 10 

NGOODsNGOOOU 

SU M Xs «!M*1X ♦ X 

X?*X* X 

SUM*? 2 

X3*x* x ? 

Siihvi *fimx3»x 3 
X4»x» X 3 

SU*"*.* *fljHyi*ML 
XLOWs A MJ MJ ( XLOW, x) 

XHI<*H= A '1 1 X 1 I X H I G M , X ) 

GO n l a 


20 I r |Nj OPO. E O. ? ) STOP "SUBROUTINE STATIST HAS NO OATA" 

N= FLO A T ( NGOOO) 

ME ANs FUK/N 
VAP=>IIMX?/N-MP AN**? 

IF(VAO.lF. 0.) Sro° "SUBROUTINE STATIST SHOWS ZE®0 VAPIANCE" 
O c V=Sf* c T ( V A A. ) 

ME AMS Q * S'J WX2 / N 
OMS=SO C T (HEANSOI 

SK c Hs (SUMX3/N- J.*MEAN*MEANSO + ?.*MEAN**3) /(OEV*VAP) 

<U®T= ( SUMX4/N- 4. *MEAN*SUMX3/N*6. *MEAN**2*SUMX2/N-3 . • MF AN**4) / 
♦ < V A ° • V A R ) 


WPITi(f,TO)TITLE,FlLE»IOPLOT,NTOTALtXMIN,XMAX,NGOOO.XLOW,XHlGH, 


♦ 



8 

9 

A 

O 

c 


ME A N, VAP, rj£ v. MEANSO. PMS. SKEW, KURT ,HST APT, HS TOP 
focmat ( im ////” RESULTS OF SUBROUTINE STATIST"//// 

" TITLE", T28.6A10//" OATA FROM FILE NO ,** ,T 24 , 1 2/ / 

" OLOT I O" * T?8 , A 10 //" TOTAL NO. OF S A MPL E S", T2 4 , 16 // 

" ALLOWA LS MINIMUM", T28.G12. 6// 

" ALLOWABLE MAXIMUM", T2fl,G12. 6// 

" SAMPLES ALLOWEO",T24,I6// 

" OTSERVEO MINIMUM", T28.G12, 6// 

" OrSEOVFP MAXIMUM", T26,G12.6// 

" mean VALUE", T24, G12. s//" V AR I A NCE" , T2 4 , G 1 2 . 6// 

" STANOARO OEVIATION", T24.G12. 6//" MEAN SOU ARE O" , T 2 8 , G 1 2 . 6/ / 
" OCCT MEAN S0UARE0",T28,G12.6//" SKE WNE SS", T2 4 , Gl 2 . 6// 

" KURTOSIS", T28.G12.6//" START T IME" , T 2 9 , A 10 / / 

" STOP TIME",T28,A10) 


RETURN 



SUn? OUT I nE MSTGRAM 


7 J/74 


OPT*l 


C ' . ’• 

OF POOR QUALITY 


PTH 4, 75060 


surojurriF mstgram if ile) 

c 

COMMON /mSTOAT/X MI N, XM At , MSECS , T ALLY « NGOOO ♦ N TO T A L 
CO-nON/PlPTOAT/XSCALE (41 • YSC ALE (4) , MS T AR T , HSTOP . T I TL P < 6) , IOPL OT 
COMMJN/ TAYT I ME /OA Y ,HPMIMSC 

TMTc-,EP PILE. T ALLY tlOOKTOTAUlOOl.PEMAINIlOO) .PASSND 
Ot MANSION FROM 1130 ), TO (100) .MESS AGE I 4) 

C 

M?CCA(;r (4) sPASSNOsO 
P£ PP» N T * 0 • 05 

MCCMSPPs T 

c 

1 BA«'O'lO* :, AS?N0* 1 

EMf-O? (?^,5,M£SSAGE) PILE.PASSNO, IOPLOT 
5 PO p “AT (••MCTG’A* 4 "til,** °UN **« 12 , IX , A 1 0 ) 

CALL P f MAO< (ME SSAGE) 

C 

00 13 T - 1 * NSEGS 
10 T ALL T ( I ) s C 
C 

RP UNO r ILc 
NGCOO*MTOT AL*C 

0ELTA=1.100001*(XMAX-XMIM)/PlOAT(NSEGS) 

c 

20 °E AO( p T L r ) x 

IP (P3P ( P I L E) .ME. 0) GO TO 40 

MTCTAL -STOTAL* l 

IP (Y. L T .xriNir.O TO 20 

TFSTs *MI *| 

c 

oo Ti r « i • msecs 

TE*T= TF ST»OELT A 
IP (X. GT . TEST) GO TO 30 
TALI* ( I) sT ALLY (I) ♦ l 
NGOOO * NGOCO* l 
GO T 3 ?0 
30 CONTINUE 
C 

’ GO n 23 
C 

40 CALL OATP(OAV) 

CALL TIME(MPMINSC) 

W®ITE (f-,P :) TITLE. PILE. 10 PLOT, NT OTA L. XM IN, X MAX, NGOOO, NSEGS. 
f OAY.HBMINSC 

5C POPM4T ctMi //" RESULTS CP SUBROUTINE MSTGRA M**/// 

1 " TITLE**, T25 , 6 A 1 0/ / “ OATA PPOM FILE NO. ** , T 25 , I 2/ / 

2 ** °LOT I C** . T ? ? . A 10 //** TOTAL NO. OF S A MPl ES** . T 2 « , 15 // 

3 ** ALLOWABLE MINIMUM”, T2fl,G12. 6// 

4 ” ALL OW A OLE MA X I MUM” , T 2* , G 1 2 . 6 // 

5 ” SAMPLES ALLOWEO”,T2fl,I5// 

6 ” NO. OP SEGMENTS”, T2«, 13// 

7 *’ BATE/TIME”, T 2 8 . A 1 0 , A 9// ) 


C 


TOTAL ( 1 ) =T ALLY (1) 

REMAIN 1 1 1 = NGOO 0- T A LL Y ( 1) 
F° CM( 1 ) = X= XM I N 


71 - 


SUBROUTINE MSTG°AM 


71/74 OPT*l 


PTN 4. 2 * 75 060 


no r o i*’.NS'r.s 

TOTAL I I) s TOTAL <1-1 ) *1411 V C II 
» p HAt m II * NGOOO- TOTAL I I) 

. 60 FROM! I I sTO I 1-1 I sX*X»OELT A 

c 

to <*r,rr.Gi »xmax 

c 

MO IT? If ,Tf| (I, FPOHCI) , TO I II .TALLY (I) .TOTAL I 1 1 , RE HA IN <1 I . I«l , NSEGS 
7 C F0OM47r* SEG. NO FOOM" , T ?9 , "TO” . T 3 0 , "S A HPL ES TOTAL". 

♦ T«;7,»RFMAINlNG"//tCCC2H , I 3. 4X.2IG12.6 .2X1 , 3 II 6,4X1 /) ) 

C 

CALL PIMTST 

c 

or TIJ’N 

C 

Eho 


ORIGINAL IZ 

OF POOR QUALITY 


kj- 

*■ , 

t 


( 




i 


t 


f . 




( 
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SII«»«0'JTINE PIHIST m/ 74 O p T« 1 ORIGINAL PAGE [3 PTN «,.?»75C60 

OF POOR QUALITY 

SUP°3UTT-|f PLMIST 
C 

C0MM0H/HtST0AT/HrSTHlN,HTS7HAX,NSEGS,T ALLY, NHG009, NMTOTAL 
CO M MON/o L rTOAT/xSCALE (u) , YSC ale (4 » ,M START, MSTOP. title C 61 , notOT 
COHHDS/TAYTIHC/OAY , HP m In SC 
INTPGPP T ALLY (100) 

OlMfJSIQN T 0 3 L IN=( U) 

OATA T r° •_ I HE/" HI *41 mum HA Y I HUH m , m VAlTO OP ,, .”S , V 

C 

CALL PLTTC U.UY f 0. .-3) 

C 

XPCA'.f ( l 1 = “ 1ST HI N 
XSCAlr (?)sMISTMAX 

CALL SCALf I XSC ALE, 10 . *2, 1) 

CALL A * l~ H . . f . , "SCALAR VALUE*** -12 .10. ,270 . * XSC A L r ( 3 ) » YSCALE ( 4 ) ) 
C 

YM AYS 0 . 

C 

DO n 1 = 1 * NS r GS 

T? PT = P L P A T ( T A L L Y ( I ) 1 
iO YH AYS AHA *1 ( YMA X , T* ST) 

C 

YSCA.P ( 1 ) r C . 

YSCALt" (2) =YMAY 

CALL S r ALc (YSCALE. 8. .2,1) 

CALL A»'30<0.,G .."SAMPLES PER I NTEPVA L” , 20 ,*., 0 0 V SC ALEC 4 I ) 

« C 

XsIHISTHIN-XSCALEIJJI/XSCALPI**) 

CALL PLTMO. ,-t.*X*3) 

c 

OELCAYtAI' (HISTMAX-HISTHIN) /PLOAT(NSEGS) 

. c 

cp ?o i = i.Hsrc-s 

YrPLDAT(tA L LY(I))/YSCALE(U) 

CALL PI OT ( Y,-l ,*X, 2) 

Y= YOfL T AX/XSC ALE (4) 

20 CALL °LOT (Y ,-l ,*X, 2) 

C 

‘ CALL PLOT (0,-1 . • X . 2 > 

CALL SYMBOL 10. . 1. 5. .15. TITLE. 0., 60 
CALL SYMBOL ( 9 . ,-5. 87 , •12,T0° LINE *270. « 31) 

CALL NMHCFPM.75,-5.97,. 12,HISTMIN,2 7 0., 3) 

CALL MUHC£ o(3. 75, -7. 35, .12, HIST HAY, 270. ,3) 

HG^OlsPLOAT (NHGOOO) 

CALL Ml'HPe o (h . 75, -A. 50 ,. 12, H GO 00 ,270., -1) 

CALL SYN1CL ( 3. 75,-9. 35,. I 2, •*/**, 2 70 . , 1 ) 

HTOTAL=PLOAT(NHTOTAL) 

CALL N(JM9rP(3. 75,-9. 50 ,. 12, HTOTAL ,270 . ,-l) 

CALL SYMBOL (3.0, -.5, . 12, I DPLOT , 270 . , 10 ) 

CALL SYHOOL(3.0,-2.,.12,MSTAor,270.,8) 

CALL S v h 10L (3, 0,-3 .12, .12,** TO**, 270 . ,2) 

CALL SYM1OL(3.0,-3.50,.12,HSTOP.270. ,8) 

CALL SYH90L (3. 0 , -7 . 75 , . 1 2 , 0 A Y, 27 0 . , 1 0 ) 

CALL SYH90L (3. G, -9 . 0 0 . . 1 2 , HR HI NSC , 27 0 . ,9) 

C 

pe T'l? N 



OneOJTlNE PLNC9H 7 1/74 OPT*l F TN 4. 2* 75060 

ORIGINAL PAGE f9 
OF POOR QUALITY 

, SUPROUTInE PLNOPH 

COMMON/STATOAT/STATMIN.STATNAX .MEAN, VAR.OEV.MEANSO.RMS, 

* S<EW,KURT,XLOW,XHlGH, NS GOOO , NST OT AL 

’ CO*MON/MlSTOAT/MlSTMlN,HlSTHAX,NSEGSf T ALLY , NMGOO 0, NK TO TAL 

CUHMON/PLCTOAT/XSCALE 14) , YSCALf !4 ) , HST ART , HST OP , T I TlE ( 6) .IOPLOT 
9£AL Hf 4N,MEANSO,KURT 
INT*‘,EP TALLY ( 100) 

PI Mf MS I ON ABSCISS! 302) .ORONATE ( 302) . TOPLINE 16) 

C 

O-TA T O p L I NE/** ME AN SOUAR~, M F RMS M . M VARIANCE'*,*' MEAN VA**, 

♦ "LUE S TO ** , ** OEV/ 

C 

x* xs: u • 1 3) 

xlNC*APSIXSCALE(4) /29.9) 

CO^P-sFLOAT (NS GOOO )* (XHIGH*X LOW) /FLOAT (NSE GS ) / ( 2 . 50 66* OE V) 

no 13 I« l » 300 
49*01 Sf ! I ) *X 

EXPON*-. 5M UX-MEAN) /OEV) **2) 

OPONA T£ ( I ) =0. 

IF<r<PON.GE.-fc75.. ANO . EX PON. LE . 7 41 . ) O RON ATE ( I) = CrSFF*EXP <EXt>ON) 

1C Xs>H I*-C 

C 

A ncrjSf ( 30 1 ) > X SC AL E ( 3 ) 

OPPNATf I lCl)«YSCALf (3) 

AP c Ct *S ( 302) s-XSCALE (4) 

O’fNA If | 3C2I *T SCALE14) 

c 

PAIL L IN" I 0R3N AT- , ABSCISS, 30 0, 1 , 0,0) 

CALL SYMBOL 19. ,1., .12, TOPLINE, 2 70. ,5 4) 

CALL NllHr»M. 75 , * 0.60 ,. 12, MEANS') ,270 . , 3) 

CALL NUMBER! A. 75,-0.70 ,. 12, RNS ,270 ., 3> 

<*ALL NUMBER! 8. 75 ,-1.85,. 1 2, VAP , 270 ., 3) 

''ALL NUMBER!*. 75,-3.25,. 12, MEAN ,270 . , 3) 
pall MIIMOfRl*. 75 ,-4. 75 ,. 12 , OEV , 2 70 . , 3) 

r 

• »F Til? N 
C 

» NO 


t 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

\ - 

10 INPUT VENTER THE NUMBER OF SCORES’* IN 
20 LET Sl«0 
30 LET S2*0 

40 PRINT "ENTER EACH SCORE (ONE AT A TIME)" 

50 FOR 1-1 TO N 

£0 INPUT X 

70 LET S1*S1«-X 

£0 LET S2«S2*X*2 

90 NEXT I 

100 LET M-Sl/N 

110 LET V»rS2-(Sl)+2/N)/N 

120 LET D-SOR(V) 

130 PRINT ,, MEAN*“MI "VARIANCE*"Vl "ST. DEV. «"D 
140 PRINT 

150 INPUT "MORE DATA (TYPE YES OR NO>"?A* 

160 IF AS=“YES" THEN 10 
170 END 


STATl : Descriptive Statistics Ungrouped Data 


10 DIM F(16> ,X<16),C<16) ORIGINAL PA: Z L7 

20 INPUT "ENTER THE NUMBER OF CLASS INTERVALS" IN OF POOR OU;' ( fTY 

7 .ET Sl-0 
4w LET S2-0 
50 LET C'Q)-0 

60 PRINT "ENTER ft FREOUNCIES. AND MIDPOINT FOR EACH CLASS" 

70 PRINT " ONE PAIR AT A TIME 

80 FOR 1*1 TO N 
90 INPUT F(I),X(I) 

100 LET S 1 *S1 ♦FC I ) *X ( I ) 

110 LET S2-S2 tF([)*X(I)t2 

120 LET Ca)»C( I-1)+F(I) 

130 NEXT I 

1A0 LET M*S1/C(N) 

150 LET V»(S2-(S1)*2/C<N))/C(N) 

1E0 LET D»SQR(V> 

170 FOR 1*1 TO N 

1E0 IF C(I)>C<N>/2 THEN 200 

190 NEXT I 

200 LET W*X(2)-X(1) 

210 LET L=X ( I ) -W/2 

220 LET Ml=L-W*<C( I-l)-C(N)/2)/F<I) 

230 PRINT "MEAN". "MEDIAN", "VARIANCE", "ST. DEV." 

240 PRINT M. M| , V. D 
250 PRINT 

2E0 INPUT "MORE DATA (TYPE YES OR NO)?"? A* 

270 IF A**"YES" THEN 20 
? END 

l 


STAT2 : Descriptive Statistics Grouped Data 


* 


QR»v*«V^-' 

OF POOR Qi 


m 


10 PRINT 

20 INPUT "ENTER FIRST SAMPLE MEAN" INI 
30 INPUT "ENTER SECOND SAMPLE MEAN" ? M2 
40 INPUT "ENTER FIRST VARIANCE" ?V1 
30 INPUT "ENTER SECOND VARIANCE" ?V2 
60 INPUT "ENTER FIRST SAMPLE SIZE" INI 
70 INPUT "ENTER SECOND SAMPLE SIZE"?N2 
80 IF NJ>«30 THEN 130 

90 LET D*SORU <Nl > *V1 -m N2>*V2> / (NI+N2-2) > *SGR( 1 /N1 *1 /N2) 

100 LET T*< M1-M2) /D 

110 PRINT "OBSERVED t«"T 

120 GOTO 170 

130 IF N2> *30 THEN 150 

140 GOTO 90 

150 LET Z=(M1-M2)/SQR<V1/N1+V2/N2) 

1E0 PRINT "OBSERVED Z="Z 

170 INPUT "MORE DATACTYPE YES OR NO)?"? A* 

IS0 'IF A$*" YES" THEN 10 
190 END 


( 


STAT3 : Testing For Significant Differences Between Two Sample Means 


ORIGINAL PAGE IS 
OF POOR QUALITY 


10 DIM M<30>, F(30> ,E(30) 

20 INPUT "ENTER THE NUMBER OF OBSERVAT IONS" 1 N 
3 INPUT "ENTER PROB. SUCC. , OR A 0 IF PROB. IS UNKNOWN"5F* 

/>RINT "ENTER ^OBSERVATIONS, #FREQUNCIES(ONE PAIR PER LINE)" 
50 FOR 1 = 1 TO N 
60 INPUT ml'.FC I) 

70 NEXT I 

£0 IF P=0 THEN 160 

90 FOR I«1 TO N 

100 LET E<I)»rtU)*P 

l 10 NEXT I 

120 LET S7*0 

170 FOP 1=1 TO N 

140 LET ST«S3-m F« I )-E« I ) >+I/E( I ) 

150 NEXT I 

1E0 PRINT "OBSERVED CHI-SOUARE VALUE *"S3 

170 GOTO 260 

1S0 LET Sl=0 

1*90 LET S2»0 

>00 FOR 1=1 TO N 

210 LET Sl*Sl+F» I) 

2.;0 LET S2-S2+MU) 

270 NEXT I 
240 LET P=S1/S2 
250 GOTO 90 

2G0 INPUT "nORE DATA (TYPE >-ES oR NOJ'V'sA* 

2/0 IF Af=" YES" THEN 20 
I'.'" END 


STAT4 : Chi-Square Distribution 


APPENDIX F 


STATISTICAL ANALYSES AND CORRESPONDING 
FALSE ALARM PROBABILITY NOMOGRAPHS 
FOR AFT I FI 6 FLIGHT DATA 


ANiSOUARE i. ..RMS ..VARIANCE CHI SQ j MEAN VALUE .! STD DEV J MI N I MUM ..MAXIMUM L .VALID DBS 



a - 1 j j t> iff rtw 






























ORIGINAL PAGE IS 
OF POOR QUALITY 

AFTIF16 FLTIA 



L-S 

L-R 

S-R 

1 

.023438 

.035156 

.011719 

2 

.023430 

-.10547 

-.12891 

3 

.0058594 

.017578 

.011719 

4 

-.017578 

.0058994 

.023438 

5 

.023438 

. 10547 

.082031 

6 

.011719 

.011719 

0 

7 

.087891 

.023439 

-.064453 

£ 

0 

.029297 

.029297 

9 

.099609 

. 13477 

.0351156 

10 

.0058594 

.041016 

.035156 

1 1 

.082031 

.029297 

-.052734 

12 

. 10547 

.029297 

-.076172 

13 

.12891 

.035156 

-.09375 

\ 

.052734 

. 10547 

.052734 

15 

-.0058594 

.22852 

.23438 

1G 

. 099G09 

. 12305 

.023438 

17 

.017578 

.041016 

.023438 

18 

-.087891 

-.076172 

.011719 

19 

.011719 

.017578 

.0058594 

20 

.029297 

.035156 

.0058594 


/J r*c 

i* 


Sensor Redundancy Management 
Report of Probability Density Function. C(DF) 
Reliability R, Probalility of false Alarm P(AF/DF) 
For the Sensor Value in the Domain of Failure 


UmI 




Parameter * 

L-S * .023438 


FLIGHT * AFT I F 1 6 FLT1A 
NORMAL ACCELERATION Sensor diffrences are L-S L-R and S-R 




L-R * .03515G 


S-R 


.011719 


P(AF/DF) Table for Sensor Diffrence .* .023438 
PDF RELIA6 . P ( AF/DF > 


ORIGINAL PAGE 13 
OF POOR QUALITY 


L-S 


. 331 394GG 


.9999 

.99969832 

. 33 1 39466 


.3983 

.99668803 

.33139466 


.3379 

.99368989 

.33139166 


.9969 

.99070374 

.33139 66 


.9959 

.98772956 

P( AF/DF) 

TABLE 

FOR SENSOR 

DIFFRENCE L-R * .03515 

PDF 


RELIAB. 

P( AF/DF) 

. 33131613 


.9999 

.99969824 

.33131613 


.9989 

.99668728 

.33131613 


.9979 

.99333838 

.33131613 


.9969 

.99070157 

.33131613 


.9959 

.98772669 

P< AF/DF) 

TABLE 

FOR SENSOR 

DIFFRENCE S-R * .01171 

PDF 


RELIAB. 

P ( AF/DF ) 

.33144173 


.9999 

.99969834 

.33144179 


.9989 

.93668849 

.33144179 


.9973 

.99369079 

.33144179 


.9969 

.99070503 

.33144179 


.9959 

.98773127 

3438 L-R 

* 

10547 S 

i-R = -.12891 

P(AF/DF) 

Table 

for Sensor 

Diffrence .= .023438 

PDF 


RELIAB. 

P< AF/DF) 

.33139466 


.9999 

.99969832 

.33139466 


.9989 

.99668803 

.33139466 


.9979 

.99368989 

.33133466 


.9969 

.99070374 

.33139466 


.9959 

.98772956 

P< AF/DF) 

TABLE 

FOR SENSOR 

DIFFRENCE L-R = -.1054' 

PDF 


RELIAB. 

P< AF/DF) 

. 33019737 


.9999 

99969722 

.33018737 


.9989 

• J9667599 

.33018737 


.9979 

.99366697 

.33018737 


.9969 

.99067007 

.33019737 


.9959 

.98768524 

P ( A^pF) 

TABLE 

liEiSI? 80 * 

DIFFRENCE p? - P/Bf7 .1289: 


( 


: 15131 


ill 

3295G19 

.3235619 

.3235619 


3373 

9363 

9359 


if 

>5 


mm 

99363303 

3303525 

9876F213 


L-3 = 


,0058534 


L-R = .017578 


S-R = .011719 


P(AF/DF) Table for Sei^so.r Diffrence .= .0058534 




fl 


RELIAB, 


P ( AF/DF ) 


.33145359 

.33145359 

.33145359 

.33145359 

.33145359 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969835 

.99668864 

.93369099 

.99070537 

.98773171 


ORIGINAL PAGE 99 
OF POOR QUALITY 


P(A^^F) TABLE ggf^gfNSOR D I FFRENCE^-^ « p j 0 1 7578 


•mm 


,33142217 

,33142217 

33142217 


9979 

9969 

9959 


liiiiii 4 

99369038 

99070448 

98773055 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE S-R = .011719 
PDF RELIAB. P(AF/DF> 


.33144179 

.33144179 

.33144179 

.33144179 

.33144179 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969834 

.99668849 

.99369079 

.99070503 

.98773127 


L-R * .0058994 


S-R * .023438 


P( AF/DF) Table for Sensor Diffrence .= -.017578 


RELIAB, 


P(AF/DF) 


.33142217 

.33142217 

.33142217 

.33142217 

.33142217 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969834 
.9966883 
. 99369036 
.99070448 
.98773055 


P < AF/DF ) TAELE FOR SENSOR DIFFRENCE L-R = .0058994 
PDF RELIAB. P(AF/DF) 


,33145352 

,33145352 

,33145352 

: 33113335 


9999 

9989 

9979 


99969835 

99668864 

99369099 




IF) TABLE 


Sensor diffrence p ^ /5fJ 


023438 


33139466 

33139466 

33139466 

33139466 

,23139466 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969832 

.99668803 

.99368989 

.99070374 

.98772956 


,023438 


L-R = .10547 


S-R = .082031 


P ( AF/DF ) Table for Sensor Diffrence .= .023438 


RELIAB 


P< AF/DF) 


,33139466 

,33139466 

,33139466 

,33139466 

,33139466 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969332 

.99668803 

.99368989 

.99070374 

.98772956 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R = .10547 
PDF RELIAE. P ( AF/DF ) 


,33018737 

,33018737 

,33018737 


.9999 
. 9969 
.9979 


,99969722 


. 9966' 


.99366697 


33018737 

33010737 


9969 

9959 


99067007 

98768524 


P( AF/DF ) 

TABLE 

FOR SENSOR 

DIFFRENCE S-R = .082031 

PDF 


RELXAS. 

P ( AF/DF ) 

33068858 


.9999 

.99969768 

33068858 


.9989 

.99668101 

33068858 


.9979 

.9936765 

33068858 


.9969 

.99068409 

33068858 


.9959 

.98770365 
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P ( AF/DF ) 

Table 

for Sensor 

Diffrence .* .011719 

PDF 


RELIAB. 

P ( AF/DF ) 

33144179 


.9999 

.99969834 

33144179 


.9989 

.99668849 

33144179 


.9979 

.99369079 

33144179 


.9969 

.99070503 

33144179 


.9959 

.98773127 


PAGE 

Quau 


P ( AF/DF ) 

TABLE 

FOR SENSOR 

DIFFRENCE L-R * .011719 

PDF 


RELIAB. 

P ( AF/DF ) 

. 33144179 


.9999 

.99969834 

.33144179 


.9989 

.99668849 

.33144179 


.9979 

.99369079 

.33144179 

.33144179 


.9969 

.9959 


P(A^F) 

TABLE 

^E.sI'.' S0R 

DIFFPENCE p S-g /5F 0 

. 3314575 


.9999 

.99969835 

.3314575 


.9989 

.99668866 

.3314575 


.9979 

.99369107 

.3314575 


.9969 

.99070547 

.3314575 


.9959 

.98773183 

087891 L-R 

= . o; 

23438 

S-R = -.064453 

P ( AF/DF ) 

Table 

for Sensor 

Diffrence .= .087891 

PDF 


RELIAB. 

P ( AF/DF ) 

.33057498 


.9999 

. 9S969754 

.33057498 


.9989 

.99667985 

.33057498 


.9979 

.99367433 

. 33057498 


.9969 

.99068092 

.33057498 


.9959 

.98769949 

P< AF/DF) 

TABLE 

FOR SENSOR 

DIFFRENCE L-R = .023438 

PDF 


RELIAB. 

P ( AF/DF ) 

.33139466 


.9999 

.99969832 

.33139466 


.9989 

.99668803 

.33139466 


.9979 

.99368989 

. 33139466 


.9969 

.99070374 

.33139466 


.9959 

.98772956 

P( AF/DF) 

TABLE 

FOR SENSOR 

DIFFRENCE S-R = -.064453 

PDF 


RELIAB. 

P( AF/DF) 

.33098259 


.9999 

.99969791 

.33098259 


.9989 

.99668393 

.33098259 


.9979 

.99368206 

.33098259 


.9969 

.39069229 

.33098253 


.9959 

.38771443 


L-R = .029297 S-R = .029297 


P( AF/DF) Table For Sensor DiFFrence .* 0 


PDF 


RELIAB, 


P< AF/DF) 


,3314575 

,3314575 

,3314575 

,3314575 

,3314575 


.9999 

.9999 

.9979 

.9969 

.9959 


. 99969935 • • 

.99663866 

.99369107 POOH QUALlTi 

.99070547 

.98773183 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R = .029297 
PDF RELIAB. P ( AF/DF ) 


,33135934 

,33135934 

,33135934 

33135934 

33135934 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969826 

.99668767 

.9936892 

.99070276 

.98772824 


P( AF/DF) TABLE FOR SENSOR DIFFRENCE S-R = .029297 
PDF RELIAB. P'AF/DF> 


.33135934 

.33135934 

.33135934 


9999 


If; 


9969 

9959 


99969826 
7 


Ilf 


99070276 

98772924 


L-S 


.099609 


L-R * .13477 


S-R * .0351156 


P( AF/DF) Table For Sensor DiFFrence .* .099609 
PDF RELIAB. P ( AF/DF ) 


3303244 

3303244 

3303244 

3303244 

3303244 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969733 
. 99667"’ J6 
.99366954 
.99067389 
.98769029 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R * .13477 

PDF RELIAB . P( AF/DF) 


32938619 

! 329386 1 9 
,32938619 


.9 999 
.9989 
. 99 7 9 
.9969 
.9959 


.99969649 


.99064757 

.98765566 


P< AF/DF) TABLE FOR SENSOR DIFFRENCE S-R = .0351156 
PDF RELIAB. P ( AF/DF > 


,33131646 

,33131646 

,33131646 

,33131646 

33131646 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969824 

.99668728 

.99368838 

.99070153 

.98772669 


L-S = 


.0058594 


L-R = .041016 


S-R = .03515G 


( 


P(AF/DF) Table For Sensor DiFFrence .= .0058594 
PDF RELIAB. P ( AF/DF ) 


.33145359 

.33145359 

.33145359 

.33145359 

.33145359 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969835 

.99668864 

.99369099 

.99070537 

.98773171 


& 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R = .041016 

PDF RELIAB. P ( AF/DF ) -87- 


.33126511 


.9999 


.99969817 


TV •' 


.33126511 
.3312651 1 
.33126511 
.33126511 


.9983 

:liil 

.9959 


.99668675 

.99772481 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE S-R * .035156 
PDF RELIAB. P ( AF/DF ) 


.33131613 

.33131613 

.33131613 


9999 

9969 

9959 


.99969824 

ripiisi 

.99070157 

.98772669 


,082031 


L-R * .029297 


S-R = -.052734 


P( AF/DF) Table for Sensor Diffrence . * .082031 


RELIAB. 


P< AF/DF) 


.33068858 

.33068858 

.33063858 

.33068858 

.33068858 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969768 
.99668101 
.9936765 
. 99038409 
.98770365 


ORIGINAL PAGE a 
0F POOR QUALITY 


P( AF/DF) TABLE FOR SENSOR DIFFRENCE L-R = .029297 
PDF RELIAB. P ( AF/DF ) 


33135934 

33135934 

33135934 


9999 


99969826 


mmn 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE S-R = -.>52734 
PDF RELIAB. P ( AF/DF ) 


3311 3952 
,33113952 
331 13952 
.33 11 3952 
33113952 


.9999 
.9989 
.9979 
.996 3 
.9959 


.99969807 

.99668549 

.99368504 

.99069664 

.98772023 


1054 7 


L-R = .029297 


S-R = -.076172 


P(A F/DF) Table for Sensor Diffrence .= .10547 


RELIAB. 


P( AF/DF) 


.33018737 

.33018737 

.33018737 

.33018737 

.33018737 


.9999 
.9989 
.9979 
. 9969 
.9959 


.99969722 

.99667599 

.99366697 

.99067007 

.98768524 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R = .029297 
PDF RELIAB. P(AF/DF> 


.33135934 

.33135934 

.33135934 

.33135934 

.33135934 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969826 

.99668767 

.9936892 

.99070276 

.98772824 


P ( AF^pF ) TABLE 


IffillS? 10 " DIFFRENCE P?fl?/DFj 


= -.076172 


33079442 

33079442 

33079442 

man 


9999 

9989 

9979 


99969778 

99668204 

9936785 

119988 * 


L-S = .12891 L-R * .033136 S-R » -.09373 

P(AF/DF) Table for Sensor Diffrence . * .12891 


ORIGINAL PAGE Q 

OF POOR QUALITY 


PDF 

RELIAE. 

P< AF/DF) 

32-5619 

.9999 

.99969662 

3295619 

.9989 

.99666969 

3295619 

3295619 

.9979 

.99365503 

.9969 

.9906525 

3295619 

.9959 

.98766213 

P(AF/gF) 

rflELE KSCiSS?* 0 * 

DIFFRENCE L-R * .035153 
P< AF/DF) 

33131613 

.9999 

.99969824 

33131613 

.9989 

.99668728 

33131613 

.9979 

.99368838 

33131613 

.9969 

.99070157 

33131613 

.9959 

.98772669 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE S-R = -.09373 
PDF RELIAe. P ( AF/DF ) 


33045358 


338 r- 


33045358 

33045358 


9999 

9969 

9953 


99969744 


205 
99067751 
98769501 


L-S 


.052734 


L-R = .10347 


S-R * .052734 


P(AF/DF) Table for Sensor Diffrence .052734 
PDF RELIAe. P ( AF/DF ) 


,33113952 

,33113952 

,33113952 

,33113952 

,33113952 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969807 

.99668549 

.99368504 

.99069664 

.98772023 


P ( AF/DF ) TABLE FOR SENSOR DIFFRENCE L-R = .10547 
PDF RELIAB. P( AF/DF) 


33018737 

33018737 

33018737 

33018737 

33018737 


9999 

9989 

9979 


.99969722 

.99667599 

.99366697 

.99067007 

.98768524 


P ( AF/DF ) TAELE FOR SENSOR DIFFRENCE S-R = .052734 
PDF RELIAB. P ( AF/DF ) 


,33113952 

,33113952 

,33113952 

,33113952 

,33113952 


.9999 

.9989 

.9979 

.9969 

.9959 


.99969807 

.99668549 

.99368504 

.99069664 

.98772023 


L-S = -.0058594 


L-R = 


,22352 


S-R = . 23438 
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NORMAL ACCELERATlbP 1 are . L -S L— R. 


ORIGINAL PAGE !S 
OF POOR QUALITY 


.98749998 .98871982 .98993966 .9911595 .99237934 .99359918 .99481902 .99603886 .9972587 .99847854 
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.995 
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ORIGINAL PAGE 13 
OF POOR QUALITY 
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ORIGINAL PAGE *3 
OF POOR QUALITY 


AFTIF16 FLT1A 


) r,r ) 

x ^ 


. 13672 

.54688 

.41016 

0 

0 

0 

-2.871 1 

-.13672 

2.7344 

-.34687 

0 

.34687 

1.9141 

1.6406 

-.27344 

-1.7773 

.27344 

2.0508 

-.27344 

-.41016 

-.13672 

.88359 

-.41016 

-1.0937 

1.3039 

1.7773 

.27344 

-.82031 

-.54687 

.27344 

.82031 

1.6406 

.82031 

1.0938 

-1.9141 

-3.0078 

.13672 

.41016 

.27344 

1 .3672 

1 . 5039 

. 13672 

.41016 

-.13672 

-.54687 

-4.2383 

-.27344 

3.9648 

.54687 

-.13672 

-.68359 

-.68359 

-.41016 

.27344 
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ORIGINAL FACfc* ’3 

OF POOR QUALITY 


Pf'< \tr « ROLL RATE 


l GYRO iirffiniti iri <.-5 / £>&'> 

.14047425 .22639666 .31231907 .39824148 .48416389 .5700863 .65600871 .74193112 .82765353 .91377594 .99969835 


6.72830956-04 


3.37512986-02 


: - \ % 
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9.99082326-02 


.1329667 


.16606517 


.19914364 


.26530058 


•29837S05 


.36453599 
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AFTIF1S FLT1A 



L-S 


L-R 


S-R 


1 

-.11713 

2 

-.20508 

3 

-.17578 

4 

-.08789: 

5 

-. 14848 

8 

-.05894 

7 

-.28387 

8 

-.23437 

9 

- . 49805 

*4 >• 

-.00789 

11 

-.29297 

12 

-.32228 

13 

.029297 

14 

-.32228 


.058594 

.058593 

0 

.11719 
. 058594 
.058593 
-.23437 
-.14848 

-.41018 

-.08789 

-.14848 

-.08789 

.23437 

-.11719 


.17578 

.28387 

.17578 

. 20508 

.0508 

.11719 

.029297 

.087891 

.08789 

0 

.14848 

.23437 

.20508 

.20508 
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QR5G1NP*L PAGE S3 

OF POOR QUALITY 


AFTIFIS FLT1C 



iudi % 




L-S 

L-R 

S-i 

1 

- . 54G87 

-.41018 

.3872 


-.82031 

-.54687 

.27433 

3 

-.41018 

-.27344 

. 13S7'2 

4 

.27344 

0 

-.27344 

5 

.95703 

-.13872 

-1.0937 

6 

-.88359 

-.41018 

.27344 

7 

-1 .0937 

-.41018 

.88359 

8 

-.27344 

-.13872 

. 13872 

9 

-.41018 

-.41016 

0 

10 

-.13872 

-.13872 

0 

11 

-.82031 

-.27344 

.54888 
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.54888 

13 
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.41018 

4-*T^ 
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0 

0 
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-.41018 

-.27344 
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Reli 



•or the Sensor value 


EfjliSI#, 

e m the OoHiin of Failure 


Par^eter * ROLL RATE OyRO^^^sor^frfnftPSr^^-S L-R and S-R 


P(f-4p/ 

.97948401 .98150544 .98352987 .9855483 .98756973 .98959116 .99161259 .89363402 .99565545 .39767688 .99969831 
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.22353154 
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